PATENT 
14XZ00088 



IN THE UNITED STATES PATENT AND TRADEMARK OFFICE 



In re Application of : METHOD AND SYSTEM FOR SIMULATING 

SUREDA ET AL. THE DIAMETER ENLARGEMENT OF A 

LESION OF A BLOOD VESSEL 

c 



Serial No. : Art Unit: ^ 



Filing Date Here\Mth : Examiner w; 

61 



: Date: November ^ , 2001 
INFORMATION DISCLOSURE STATEMENT 



[ X ] 1 . Pursuant to 37CFR 1 .97(b) 

[within three months of filing or prior to first Office Action] 

[ ] 2. Certification Pursuant to 37 CFR 1 .97(c) 

[before Final Office Action or Allowance] 

[ ] 3. Fee Payment Pursuant to 37 CFR 1 .97(c) 

[before Final Office Action or Allowance] 

[ ] 4. Petition, Certification & Petition Fee Payment Pursuant to 37 CFR 1 .97(d) 

[before issue fee payment] 



ATTN: BOX PATENT APPLICATION 
Assistant Commissioner for Patents 
Washington, DC 20231 



Sir: 



The following are submitted In the above application in compliance with 37 CFR 1 .97 and 
37 CFR 1.98. 



[ X ] 5. A list of documents on Form PTO-1449 together with copies of each identified 

document and a translation thereof or a concise explanation of each non-English 
language document or a Search Report or communication from a non-US patent 
office or an International Search Report from an International Searching 
Authority for a patent application filed via the Patent Cooperation Treaty or 
document (s) cited in the application or the priority application. 



This paper Is submitted in accordance with: 



[ X ] 6. 37 CFR 1 .97(b): [within three months of filing or prior to first Office Action] 



CERTIFICATE OF MAILING 



I hereby certify that this correspondence is being deposited with the United States Postal Service 
as First Class Mail in an envelope addressed to: Assistant Commissioner for Patents, 
Washington, DC 20231, on the date indicated below. 

B y IpS^*^ ^ ^ Date: November 9 . 2001 

Dolores Tillson 

Page 1 of 2 



[ ] 7. 37 CFR 1 .97(c): [before Final Office Action or Allowance, whichever is earlier]; 
and 



[ ] 8. The required Certification made in item 1 1 below; or 

[ ] 9. The $240.00 fee specified in 37 CFR 1 .17(p) for submission of this 
Information Disclosure Statement is authorized in item 14 below. 



[ ] 10. 37 CFR 1 .97(d): [before issue fee payment]; and 



a) This Is a petition for consideration of the subject Information Disclosure 

Statement. The petition fee of $130.00 required by 37 CFR 1.17(i)(1) is 
authorized in item 14 below. (Direct this letter to: "Attention PETITIONS 
EXAMINER" and if applicable include batch locator information, e.g., "Allowed 
Files, Batch , Date of Allowance and 



b) The required Certification is stated in item 1 1 below. 



[ ] 1 1 . Certification 



[ ] 12. Each item or information contained in this Statement was first cited in a Search 

Report or communication from a non-US patent office or an International Search 
Report from an International Searching Authority for a patent application filed via 
the Patent Cooperation Treaty in a counterpart foreign patent application not 
more than three months prior to the filing of this Statement; or 

[ ] 13. No item contained in this Statement was cited in a Search Report or 

communication from a non-US patent office or an International Search Report 
from an International Searching Authority for a patent application filed via the 
Patent Cooperation Treaty in a counterpart foreign patent application or, to the 
knowledge of the person signing this document after making reasonable inquiry, 
was known to any individual designated in 37 CFR 1 .56(c) more than three (3) 
months prior to the filing date of this Statement. 



[ ] 14. Payment of all applicable fees: 



D Please charge all applicable fees associated with the submittal of this 

Information Disclosure Statement to Deposit Account No. 09-0470. An original 
and two (2) copies of this document are enclosed. 

D Enclosed is a check in the amount of $ in payment of all applicable 

fees associated with the submittal of this Information Disclosure Statement. An 
original and two (2) copies of this document are enclosed. 



{esppctfully submitted, 




Address: General Electric Company 

3135 Easton Turnpike, Fairfield. CT 06431-0001 
Fax: 203-373-3991 
Tel: 203-373-2867 

Page 2 of 2 



L h 









MPOAMATION DiSClOtUM CItATlOM 






9^hm% MM 






■5 mwm%f 



OtMA MCWMlWTt (Hm^,m4 Ammm. r*i*. ^« 

























pMLridrLU/ odchvD^o <m^\ i{rftsfi^aof>tp(mH<^ /N^rn/rc^j f^t^e^ ffuui^^^ 

OtJUu^ ^OL AlO. 'J SJ^ _ 








^^^^K^^^ ^2;:sr P.)ACy^r^. g^/^r^P /nai^mr//^i4£D 







•CBA«iM I. teMi It cmiffi tMi^^W. 



(HIV. 

INFORMATION DiSCiOSUM CITATION 
(Vm m^mmt t^9^t% if m9m 9 m^t 

at. MTI 








9^tm% MM 





19 *»••••• 



DEMANDE DE : BREVET 

NO : 0015691000 DU 04/12/00 

V/REF. : B99/2540FR AK/G 




INSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



CASALONGA 



BUREAU CASALONGA ET JOSSE 

8 AVENUE PERCIER 
75008 PARIS 



PARIS, LE 16 AOUT 



2001 



OBJET : NOTIFICATION D'UN RAPPORT DE RECHERCHE PRELIMINAIRE 
AVEC REPONSE OBLIGATOIRE 



Messieurs, 

J'ai I'honneur de vous adresser, en annexe, le rapport de recherche preliminaire etabli conformement a 
I'articie R.612-57 du code de la propriete intellectuelle, citant les documents qui peuvent etre pris en consi- 
deration pour apprecier la nouveaute et I'activite inventive de {'invention, objet de votre demande. 

Salon Tarticle R. 612-59 du code precite, vous disposez d'un delai de 3 mois a compter de la date de 
reception de ce rapport de recherche preliminaire pour y repondre par ecrit. Avant Texpiration de ce delai, 
celui-ci peut etre renouvel6 une fois sur votre requete. 

Suivant la categorie des documents cites, vous pouvez etre tenu a une obligation de reponse (par exemple, 
si le rapport de recherche preliminaire mentionne des documents de categorie X ou Y). Dans ce cas, un 
papillon rouge est appose sur cette lettre et le defaut de reponse entramera le rejet de la demande. Dans 
le cas contraire, ce papillon est jaune. 

Dans tous les cas, il est de votre interet en elaborant votre reponse, de tenir compte de tous les 
documents cites. 

Selon les articles R. 61 2-58 et R. 61 2-60 du code precite, votre reponse peut consister : 

- soit en de nouvelles revendicatlons (en 3 exemplaires). Dans ce cas, vous devez signaler les 
changements apportes aux revendicatlons initiales. Vous pouvez y joindre des observations qui mettent en 
evidence les caracteristiques techniques de ces nouvelles revendications qui echappent a Topposabilite 
des anteriorites citees, 

- soit seulement en des observations qui ont alors pour objet de discuter Topposabilite des anteriorites 

citees. 

Veuillez agreer I'expression de ma consideration distinguee. 



Pour le Directeur general de I'Institut national 
de la propridt^ industrielle 



Le Chefdu d^partement des brevets 




SIEGE 

INSTITUT Petersbourg 
OF 75800 PARIS cedex 08 
NATIONAL ot UlAphone : 33 (1) 53 04 53 04 
LA PROPRJETE T6ttcopte : 33 {1M2 93 59 30 
I N 0 U S T R r e L L E www.inpi.fr 



55S/090501 



ETABLISSEMENT PUBLIC NATIONAL CREE PAR LA LOI N* 51-444 DU 19 AVRIL 1951 



IN'I 



IINSTITUT 
NATIONAL DE 
LA PROPRIETE 
INDUSTRIELLE 



RAPPORT DE RECHERCHE 

PRELIMINAIRE 

etabli sur la base des dernieres revendications 
deposees avant le commencement de la recherche 



N" d'enreglstrement 
national 



FA 597502 
FR 0015691 



DOCUMENTS CONSIDERES COMME PERTINENTS 



Cat^goiie 



Revendication(s) 
concernee(s} 



Citation du document avec indication, en cas de besoin. 

des parties pertinentes 



BALAKRISHNA HARIDAS, CLINTON A. HAYNES: 

"Predictive Analysis at the Forefront of 

Medical Product Development" 

MEDICAL DEVICE & DIAGNOSTIC INDUSTRY 

MAGAZINE, 'en ligne! 

octobre 1999 (1999-10), XP002173874 

Extrait de 1 'Internet: 

<URL : http : //devi eel i nk . com/mddi /archi ve/99 
/10/009.html> 'extrait le 2001-07-30! 

* le document en entier * 

"Simulation of the In V1vo Deployment of 
a Palmaz-Schatz Type Stent 1n a Partially 
Occluded Artery" 

SUPERC0MPUTIN6 INSTITUTE, RESEARCH 

BULLETIN ONLINE, 'en ligne! 

vol. 14, no. 3, juillet 1998 (1998-07), 

XP002173875 

Extrait de 1 'Internet: 

<URL : http : //web . ms 1 . umn . edu/gener al /Bui 1 et 
in/Vol .14-No.3/j uly98.html > 
'extrait le 2001-07-30! 

* le document en entier * 

HOLZAPFEL G A ET AL: "LARGE STRAIN 
ANALYSIS OF SOFT BIOLOGICAL MEMBRANES: 
FORMULATION AND FINITE ELEMENT ANALYSIS" 
COMPUTER METHODS IN APPLIED MECHANICS AND 
ENG I NEE RING , NORTH-HOLLAND , AMSTERDAM, NL , 
vol. 132, no. 1/02, 

15 mai 1996 (1996-05-15). pages 45-61, 

XP000972573 

ISSN: 0045-7825 

* page 56 - page 59 * 

-/-- 



1-7,9-11 



1-7,9-11 



1,11 



Classement attribu^ 
a ('Invention par I'INPI 



A61 86/00 

A61B19/00 

G06F17/00 

606T17/00 

G06T17/20 



OOMAINES TECHNIQUES 
RECHERCHte (lnt.CL.7) 



G09B 



Date cf achfevement de la recherche 

3 aoQt 2001 



Examinateur 

Szarowski, A 



CATEGORIE DES DOCUMENTS CITES 

X : parttculierement pertinent a lui seul 

Y : particulierement pertinent en combinaison avec un 

autre document de la m6me categorie 
A : arriere-plan lechnologlque 
O : divulgation non-6crite 
P : document intercalaire 



theorie ou principe d la base de {'invention 
; document de brevet b6neficiant d'une date anterieure 

^ la date de dep6t et qui n'a 6t6 publie qu'6 cette date 

de depdt ou qu'^ une date post^rleure. 
: cite dans la demande 

cit6 pour d'autres raisons 



& : membre de la m6me f ami He. document correspondant 



page 1 de 2 




INSTITUT 
NATIONAL DE 
LA PROPBICTE 
INDUfiTRIELLE 



RAPPORT DE RECHERCHE 
PRELIMINAIRE 

etabli sur la base des dernieres revendications 

deposees avant le commencement de la recherche 



N" d'enregistrement 
national 



FA 597502 
FR 0015691 



DOCUMENTS CONSIDERES COMME PERTINENTS 



Categorie 



Citation du document avec indication, en cas de besoin, 
des parties pertinentes 



Revendicat(on(s) 
concern^e(s) 



Ciassement aUrtbu^ 
a I'invention par I'INPI 



D,A 



6. FINET: "Echographle endocoronaire et 
angioplastie" 

ARCHIVES DES MALADIES DU COEUR ET DES 
VAISSEAUX, 

vol. 92, no. 11, novembre 1999 (1999-11), 
pages 1681-1689, XP001014246 
* le document en entier * 



1,11 



DOMAINES TECHNIQUES 
RECHERCH^S (lnt.CL.7) 



Date d'achevement de la recherche 

3 aoQt 2001 



Examinaleur 

Szarowski, A 



CATEGORIE DES DOCUMENTS CITES 



X : paniculterement pertinent a lui seul 

Y : particuiieremenl pertinent en combinaison avec un 

autre document de ta m6rne categorie 
A : arrifere-plan technoiogique 
O : divutgation non-^crite 
P : document intercalaire 



theorie ou principe a la base de I'invention 
; document de brevet b6n6fidar»t d'une dale ant^rieure 

^ la date de d6p6t et qui n'a ete publie qu'd cette date 

de dep61 ou qu'^ une date posterieure. 
: cite dans la demande 

cite pour d*autres raisons 



& : membre de la m6me famille. document correspondant 



page 2 de 2 



■XP-0021 73874 



^ ^ dfreetoiies Inffpgtry |>i^sgCenter commur^^ search 



Si 



if industry 



MedfcaJ Device & D^gnostic Indust/y Magazine 
MODI Article Index 



An M7ifZ7/ October 1999 Column 



DESIGN 



P.D. OCX [^Jnfj 

P. /'/o 




Predictive Analysis at the Forefront of Medical 
Product Development 

Fredictfve ansfy^'s and srmu/ai/on ofwecf/ca/dey'fce and /mp/ant pe/fyrmanoe as/ng no/j//near 
finite element anafysis can s/gnWcandy accelerate pnx/uct devefopment cyc/es whife he/p/ng 
mana^sfctu/ers to avofdcosdym/stakes. 

Balakrishna Haridas and Clinton A. Haynes 

The biomedical device and implant industry provides a unique and challenging environment for 
medical product development. Engineers are actively involved in the design, development, and 
manufacture of materials and devices intended for a variety of life-sustaining applications. 
Because the functions of such devices and implants can have varying levels of impact on the 
v/ell-being of patients, FDA has three levels of classification for regulatory purposes — Class I, 
Class II, and Class III. The increasing levels of regulation from Classes I through til represent, in 
eii&cX, a rigorous quality control program to ensure that new devices undergo adequate and, if 
required, extensive testing and trials to establish their performance prior to formal release into 
the U.S. market. 

Understandably, the development cycle for a medical product can be quite complex, especially 
for Class 11 and Class III devices. The conventional product develop-ment process of 
design-build -test-re build typically involves several trial-and-en-or-based iterations of a concept 
prior to full-scale production. These iterations are costly and tax resources, especially since 
rework loops in the initial development phases of the device do not directly verify its in vivo 
function. Popular rapid prototyping technologies such as stereo-lithography, fused deposition 
modeling, and selective laser sintering allow design teams to duplicate part geometry and 
qualitatively verify design intent. Unfortunately, rapid prototyping does not provide any 
quantitative feedback on the functional mechanical performance of an actual device in vivo. 

In the early stages of the development of a new device or implant, it is not desirable 
economically to manufacture functional prototypes of several concepts for in vivo testing. When 
prototypes are fabricated, the inability of a device design to deliver the intended function is 
established only after the organization has gone through an entire product developn-.ent cycle. 
In this case, a change in the design could have a dramatic financial impact on the company that 
is developing the product. 

Current approaches to medical device development can be significantly enhanced using such 
advanced computational methods as nonlinear finite element analysis (FEA). This 
computer-t)ased simulation technology has experienced tremendous growth over the past 10 
years. As a result, there are several commercial sottv/are programs availat>le to the engineer 
appropriately trained in the theory and practical application of this analysis method. 

The theoretical simulation techniques descrit)€d below have a definite role in the engineering 
development cycle for medical devices or implants. While challenges do exist in the application 
of the technology, an engineer with the appropriate level of experience and training can 
effectively apply these methods to accelerate the product development process. 

Predictive analysis will never completely replace prototyping and testing. However, drawing 
from past experience with this technology, in most cases predictive analysis tools limit the 
number of functional prototypes and trial testing, including expensive animal trials.^ Effective 
integration of engineering-mechanics-based design decisions into the medical device 
development process is also usually achieved. 

THE ROLE OF PREDICTIVE ANALYSIS AND SIMULATION 

Two popular yet diametrically opposed perspectives on the intended role of predictive 
simulation in the nriedical products industry are both based on misconceptions. The first school 
of thought has reservations regarding the application of computer-based simulation. The 
reasoning behind this is that inputs to a computational simulation of the mechanics of a device 
and sun-ounding tissue sometinries appear to misrepresent the true boundary conditions or 
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material properties. The net result of this misconception is that engineering analysis and 
simulation as a tool is abandoned, and a trial-and-error process is adopted as the preferred 
product development strategy. The second school of thought believes that simulation will 
completely replace the need for all prototyping. While this has occurred to a limited extent in 
conventional engineering applications such as the automotive industry, it is not likely to occur 
within the medical device industry because of the kind of products involved. 

The true role of computational simulation and predictive analysis lies between these two 
extremes. A carefully defined computational mechanics model of a device-tissue interaction 
problem can be invaluable to a biomedical designer in providing feedback on the relative merits 
and demerits of different design options, for example. In this sense, simulation of the 
mechanical performance of a medical device or implant in vivo using nonlinear FEA can be 
employed routinely at different phases of the product development cycle as deenned appropriate 
by a project team. 

Most medical products involve materials, such as polymers, and manufacturing techniques 
(e.g., injection molding, extrusion, blow molding, thermoforming, and compression molding) that 
introduce significant variat>ility in component performance. Designs must also account for 
biological variability in cases where the device interfaces directly with tissue. In this situation, 
highly variable tissue material properties play a significant role. Mechanical properties of living 
tissue are functions of several factors such as age, pathological state, type of tissue, and 
collagen-elastin-muscle content. The resulting device-tissue interaction problems involve severe 
nonlinearities in mechanical and material behavior. This is the underlying feature of most design 
problems involving polymeric components that interact with tissue, which is a highly nonlinear 
material. As a result, some form of in vitro or in vivo testing in animal models with a fully 
functional prototype is always required. Nonlinear FEA can be used as a complementary tool In 
conjunction with appropriate material descriptions to predict such nonlinear phenomena. This 
adds tremendous value to the product devebpment process and dramatically reduces the risks 
associated with nnechanical component design. 

APPLICATIONS OF PREDICTIVE ANALYSIS AND SIMULATION 

The growth and development of the application of nonlinear FEA technology as a predictive 
analysis tool has enabled its application in several areas of medical device design and 
manufacturing? The diverse range of applications include such hard tissue implants and 
devices as: 



• Fracture wires. 

• Fixation pins (internal/external). 

• Bone plates and screws. 

• Total joint replacements. 

• Dental implants. 

Soft tissue implants and devices include: 

• Vascular grafts. 

• Percutaneous devices. 

• Heart assist devices (pacemakers and LVADs). 

• Ear and eye implants. 

• Surgical tools. 

As mentioned above, the role for predictive analysis tools can encompass several product 
development functions, such as upstream design concept evaluation, "what-if material 
sensitivity studies, understanding tolerance effects, manufacturing and materials processing 
simulation, feilure analysis in n>eta!licand polymeric components, and generation of analytical 
data in support of 510(k) and pre market approval applications. 

Applications of predictive analysis for the first three product design and development activities 
from the above list are described below. A forthcoming article will address the application of this 
technology to the design and optimization of manufacturing processes, as well as to problems in 
failure analysis in metallic and polymeric components in medical devices. 



UPSTREAM DESIGN CONCEPT EVALUATION 

One of the most effective applications of predictive simulation is upstream design concept 
evaluation in the earliest stage of the product development cycle, i.e., when a designer is 
conceptualizing the configuration of the device. Application of predictive nonlinear FEA at this 
stage facilitates physics-based quantitative ranking of several design concepts. The following 
case study is a good example of the application of FEA to these areas. 



i 




F/gure 1. Sche/n3t/c of a typ/ca/ 
percutaneoLfS de}/ice. Ax/a/and 
torsJona/ mot/ons oouJd devBk>p 
stresses in the devfce as we// as at 
the tfssue-//np/a/7t //iterfaoe. 



Case Study: Percutaneous Implant Evaluation. An example of the application of predictive 
analysis is the design evaluation of a conceptual percutaneous implant (Figure l).*^ 
Percutaneous implants are surgically introduced under the skin for long-term use.^ These 
implants serve as conduits for the transfer of information in applications that can include drug 
delivery, blood monitoring, colorectal surgical procedures, pacemakers, and ports for kidney 
dialysis. Mechanical design issues for such devices include: 

• Material selection. 

• Minimization of tissue-device interface stresses. 

• Optimization of implant geometry for various possible motions of the device. 

Figure 2. Frnite e/e/nentmode/af 
(he percutarwous devfCB, s/x>mt 
in F/gure 1, fmp/anted ^ skrn. 





: Figure 3, Von M/ses stress /n the 
\ device during torsionai k>admg 
: andniotAon. 



A key issue in designing these implants relates to the local stiffness of the implant at the 
implant-tissue interface. A large mismatch in stiffness between the implant surface and the 
tissue will result in stress concentrations that could be detrimental to the implant material, the 
tissue, or both. Given that the properties of soft or hard tissue are very sensitive to changes in 
mechanical stimulus in vivo, a phenomenon known as Wolffs Law.^the engineering of the 
tissue -implant interface is extremely important. A predictive simulation rTx>del allows the 
performance of the design to be bounded in this sense. Nonlinear FEA results for the 
percutaneous implant example are shown in Figures 2 through 5. These simulations capture 
interiaminar stresses within the implant materials in addition to obtaining estimates of the loads 
transferred to skin across the Dacron-tissue interface. 




F/gure 4. Von M/ses stress /n the device durf/tg ax/a//oad/ng and motk>n. 




F/gure 5. Stress m the sk/n Jayers during ax/a/ motfon of the deyt'ce. 

Results from the nonlinear FEA simulation model of this percutaneous device can be used to 
drive changes in the design, vs/hlch could Include material selection or geometric changes that 
would ensure the long-term reliatnlity of the device and minimize the mechanical stresses at the 
tissue-device interface. The efficiency of the seal at the port can also be evaluated with FEA 
and contact-stress analysis to determine the magnitude of the seal pressure. Because the seal 
is microbial, it is critical that a minimum contact pressure is maintained over time and across the 
range of manufacturing tolerances in the access port. Recently described methods of using FEA 
to analyze sealing and closure systems can be applied to this class of design problems. 

•WHAT-IF- MATERIAL SENSITIVITY STUDIES 

The unique feature of designing a medical device for a load-bearing application is the fact that 
the device material chosen not only has to be biocompatible and FDA approved, but also must 
deliver the appropriate mechanical function within the design envelope. A computational model 
can be used to perform "what-if evaluations to determine the optimal range of mechanical 
properties of the device material. This simulation can be leveraged to assess the effect of 
biological variability and resulting mechanical property differences in the tissue that could Impact 
device performance. The measurement of the tissue mechanical properties, which are often 
highly nonlinear, can be characterized using specialized in vitro tests on cadaveric tissue, in 
vivo tests in animal tissue, or a combination of both.^ Sample sizes for these tests can be 
controlled via statistical power analyses to deliver ±2o, ±3a, or even ±60 range of properties 
(o — standard deviation).® Such material data allow the simulation of the device-tissue 
interaction problem to be exercised over the entire range of statistical variability, allovrtng the 
function of the device and the materials to be fully assessed. 

Case Study. Balloon Angioplasty Catheter Evaluation. Percutaneous transluminal 
angioplasty (PTA) is a surgical procedure performed to alleviate stenosed arteries by 
compressing the calcified plaque that can Nock the flow of blood.^ The goals of this procedure 
are to enlarge the lumen of the stenosed vessel, to maintain the lumen over time, and to provide 
an intimal surface that promotes neo-intimal hyperplasia. 
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Figure 6. Pressure vs. d/arrieierra/atk^nsh/p for a typical ba/kx>rj catheter during fi^^ 
cond0orjs, 

A critical design variable that affects the performance of a balioon catheter in terms of geometry 
and material selection is its bifurcation pressure and subsequent inflation. The bifurcation 
pressure for an inflating membrane is the internal pressure required to initiate dilation, beyorul 
which the growth of the membrane proceeds in an uncontrolled fiashion (Figure 6). The 
mechanics of memlxane inflation fall under a particular class of highly nonlinear phenornena 
called snap-through-bJfurcation problems. Typical pressure vs. radius curves, illustrated In 
Figure 6 for balloon catheters, can be computed with a finite element model of free inflation 
(Figure 7) that incorporates the appropriate nonlinear material data for the polymer or elastomer 
used in the catheter. In general, these catheter materials can be described by hyperelasUc or 
plasticity material models that are available in most leading nonlinear finite element codesJ The 
characterization of these materials, especially hyperelastic materials, has to be done carefully 
because a standard uniaxial tension test alone will not yield adequate data for balloon inflation 
simulations. Additional tests to characterize the response of the polymer under constrained 
conditions, such as biaxial stretching, must be performed to generate the appropriate material 
constants for input into the finite element calculations. 

Prior knowledge of the bifurcation pressure is extremely useful to the designer of an angioplasty 
balloon system because the inflation of the membrane progresses in an unstable fashion 
beyond the bifurcation point. This can be seen in the form of large radial expansions with no 
increase In cavity pressure (Figure 6). The PTA balloon manufacturer would thus have useful 
information available regarding the window of operation of a given catheter and balloon design. 
Design variables, including material properties, wall thickness distribution, and geometry, could 
t>e used to control the initiation of the bifurcation and subsequent inflation characteristics of a 
given device. 
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Figure Z FEA simuiation of free inffation oftaiioon cat/feterto determine tmircatidn pressure. 
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Figure S. 3-D ffn/te e/ement mode/ofstenosecfarte/y m'£h /Sxedmre itaffoon caiheter prior to 
/nffatforj. 

Once the bifurcation pressure has been determined, the optimal balloon geometry and wall 
thickness can be determined by performing simulations of the surgical procedure. This requires 
finite element modeling of the balloon positioned within the lumen of the stenosis and then 
inflated against a virtual blocked vessel (Figures 6-10). This analysis requires characterization 
of the 3-D geometry of a stenosis/lesion using imaging and x-ray data, plaque material 
properties ranging from highly calcified to minimal calcification levels, and hyperelastic material 
properties for a normal artery through appropriate material testing. 

As a starting point, this information is typically obtained from published data and research 
articles and is usually sufficient for design evaluation purposes. Such a virtual predictive 
simulation model is an ideal platform to generate catheter performance data in terms of the 
ability of the design to compress plaque without initiating tears at the plaque-intine interface. 
Intimal splitting of plaque is clinically detrimental to the success of the angioplasty procedure.® 
Evaluation of the catheter in this virtual model can also l>e exercised over a ±6a range of 
plaque and artery properties as described earlier. The introduction of a supporting metal stent 
can also be simulated in the same model. In this case, the goal is to design in the stent the 
appropriate radial stiffness that would te required to support the dilated arterial cross section. 

UNDERSTANDING TOLERANCE EFFECTS 

Manufacturing tolerances have a significant impact on the performance of the many medical 
devices and implants that use polymeric materials. If the range of dimensional tolerances that 
are typtcal for the process is known, a predictive analysis model of the device can be assessed 
at the nominal as well as at the extremes of the expected dtmerisional variability. Analytical data 
generated from such models can be used in Monte Cario simulations to generate the statistical 
range of the performance displayed by the device or implant. 

Case Study: Evaluation of a Blood Vessel Ligation Clip System. Ligation clips are being 
used increasingly during surgical procedures as a means to temporarily ligate blood vessels that 
feed into the surgery site. These clips are usually loaded in cartridges and are dispensed Into 
position by an actuating mechanism that drives a stack of clips in the cartridge through a rigid 
metallic channel (Figure 11). As the device is fired, the clips pass through a srriall bump-like 
feature in the channel called a retention bead. The force exerted by the surgeon to move the 
clips through the device is a direct function of the interference between the clip dimension (Z^ 
and the retention bead dimension {(CT). Depending on how the clips and channel were 
manufactured, there will be tolerances associated with the resulting interference that must be 
overcome during device actuation. These tolerances directly affect the force input that must be 
exerted by the surgeon. For example, metal clips formed by a bending operation will have 
dimensional variations caused by the varying degrees of elastic spring-back following the 
forming operation. Injection-molded clips, on the other hand, will exhibit variability caused by 
postmolding shrinkage and mold-tooling-induced tolerances. 
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A nonlinear FEA model of thi& device idealized as a three-clip system can be used to 
understand the effects of clip tolerances on the firing-force profile, as well as clip-to-clip 
Interactions and stresses (Figure 12). The firing-force time ounces descritjed in Figure 13 are 
calculated by the simulations and describe the variability in the device firing force as a function 
of dimensional tolerance and of the stacking sequence of the clips (minimunvmaximum-nomirml 
in this case). 

The unique aspect of this finite element model of the clip system is that there is a strong 
coupling between inertlal and deformational effects in the clips. To fully capture this coupling, 
the model has been solved as a structural dynamics problem. A plot of the different energy 
quantities as a function of time, during firing of the three clips, deariy shows the energy 
accumulation from straining (i.e., deforming) each clip as rt passes through the retention bead 
(Figure 1 4). )t should be emphasized that almost all of this strain energy is converted into fdnetio 
energy as each clip "jumps out" and moves downstream from the retention bead. 

Another interesting result is that the force to fire the maximum-dimension clip is slightly lower 
than the value for the nominal-dimension dip (Figure 13). This is explained by the fact that the 
nominal clip slightly compresses the maximum clip because of c!ip-to-dip interaction during 
actuation. As a result, the effectK^e dimension D is reduced to a value below the nominal at the 
point when the maximum clip passes through the retention bead, thus dropping the firing force 
of the clip to a value below the nominal. 
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Results from analyzing several stacking sequences of this three-clip system can be used to 
develop statistical models of clip firing to fully understand the effects of an art)itrary number of 
clips in the device. Once the clip is positioned around a blood vessel, an elastic-plastic finite 
element simulation can be used to calculate the force needed to achieve ligation In a Wood 
vessel of given geometry and matenal properties (Figures 16 a~d). This solution can be used 
subsequently to either optimize the mechanical advantage of the finng mechanism, or to 
redesign the clip geometry to minimize the forming force and enhance the effectiveness of the 
ligation mechanism. 

TECHNOLOGICAL CHALLENGES 

Effective application of predictive analysis requires a clear understanding of the technological 
challenges associated wrth developing reasonable computational models of in vivo device 
function. The most significant issues routinely encountered include estimation of device loads 
and tx)undary conditions, proper representation of tissue geometry; proper representation of 
device and tissue material behavior, and quantitative validation of predictions with controlled 
small-scale tests. 
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Project teams considering the effective use of predictive analysis in the product development 
process should ensure that the above issues are adequately handled within the project timeline, 
issues related to estimation of device loads and boundary conditions can be addressed by 
ensuring that an appropriate level of understanding has been achieved with regards to the 
design intent of the device. Biomechanical principles can be employed to derive realistic 
fc>oundary conditions of the device or tissue for FEA model representation.^^ Tissue geometry 
and composition issues can be quantitatively addressed by appropriate imaging niodalities, 
including computer-aided tomography* positron emission tomography, and x-ray, as well as 
microscopy techniques (e.g., histology, morphometry) .'^^•'^^ Device material response, in 
general, can be handled via conventional material4esting methods, tf loads transmitted through 
tissue are important, however, characterization of tissue response will require fairly specialized 
testing in vivo or in vrtro. Statistics play an important role in the analysis of tissue material data 
because of the need to understand variability of the properties at a desired level of significance . 

CONCLUSION 

Predictive analysis or simulation of medical devices and implants is a product development tool 
that can significantly accelerate the time to market and help manufacturers avoid costly 
mistakes early in the design process. The examples discussed atx^ve demonstrate the 
successful application of nonlinear FEA for evaluation of medical devices. They incorporate 
several of the technical challenges that are commonly faced in medical implant and device 
design. Although these technological challenges are inherent to the functional and regulatory 
requirements placed on medical products, they can be handled by carefully designed 
experiments and the application of biomechanical engineering principles. 

Simulation methods offer the medical product designer an opportunity to explore the fijnctional 
attributes of a design before tooling and trial-testing commitments are made. More importantly, 
the results of these simulations provide quantitative data and insights that are not available via 
prototyping methods alone. Issues such as manufacturing defect sensitivity, dimensional 
variability, material selection, and variability associated with user operation can be addressed to 
evaluate the robustness of a particular design concept. 

With increasing demands being placed on medical device manufacturers by the market and by 
regulatory agencies, it is critical to establish a robust product development process that 
incorporates sound engineering. This can be achieved by incorporating a predictive engineering 
Simula tion-Hdriven design evaluation process within the design loop. Drawing from past 
experience, such an integration of predictive simulation into the design cycle has always 
produced tangible benefits in terms of cost reduction, reduced time to market, and minimal 
rework loops. 
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Simulation of the In Vivo Deployment of a Palmaz-Schatz 
Type Stent in a Partially Occluded Artery 

in this is sue 

Balloon angioplasty has proven to be an effective alternative to more invasive 
surgical procedures for the treatment of coronary disease, the leading cause of 
death in the United States. The goal of balloon angioplasty is to permanently increase 
the luminal size of a blood vessel that has been expanded by the inflation of a balloon at 
high pressure within the narrowed vessel segment to mechanically remove plaque 
deposition responsible for occlusion. 

During angioplasty, a mechanical device known as a stent may be deployed at the site of 
vascular injury. Stents may be fabricated from a number of biocompatible materials, 
typically metals, and are preferably radio-opaque. Optimal stent deployment technique 
leads to uniform circumferential expansion and very close apposition to the vessel wall. 

Unfortunately, balloons and stents are not designed for deployment as matching systems. Physicians often 
rely on personal preference and experience in selecting a combination of angioplasty balloons and stents. This 
may lead to inconsistent expansion of system components against the vascular wall. Stent deployment in an 
irregularly occluded or eccentric lumen can be particularly problematic. 

Professor William Gleason in the Department of Laboratory Medicine and 
Pathology at the University of fs/linnesota and Svenn Borgersen of 
BIOSIMulations Inc., Eagan, fVlinnesota used a Palmaz-Schatz type stent 
geometry as the basis for development of a three dimensional finite element 
analysis (FEA) model to simulate stent deployment within a partially occluded, 
asymmetric arterial lumen as shown in the figures. The purpose of the analysis 
was to examine the effects of an asymmetric lumen on deployed stent 
geometry and stress levels, 

undepioyed stent eccentric arteriaiiumen Modol development and analysis were basod on the mentat/marc k6.2 fea 

software from MARC Analysis Research Corporation, Palo Alto, California. The 
model consisted of eight noded, isoparametric, hexahedral elements. This element type has three translational 
degrees-of -freedom per node and a total of twenty-four degrees-of-freedom per element with eight Gaussian 
integration points. The model allows for large displacement, geometric stiffening, non-linear material 
properties, and full three dimensional contact. For this study, vascular walls and lumen occlusion surfaces 
were represented as three dimensional mathematically rigid boundaries. In order to evaluate anticipated stress 
concentration effects due to geometry at junctions of longitudinal and circumferential stent members, model 
mesh refinement was included at these locations. 

Expansion of the stent to its maximum design configuration was accomplished using a uniformly distributed 
loading by use of a follower force on the internal surfaces of the stent. The magnitude of the pressure load was 
increased incrementally until the stent conformed to the simulated vascular geometry or the ultimate strength 
of the material was exceeded. 

The obtained results clearly illustrated the inherent risks of deploying this type of stent in an asymmetric 
lumen: asymmetric deployed geometry, crippling of the structure, and non-conformity of the stent to the 
vascular lumen. A stent deployed under these conditions could be succeptible to collapse. Lack of structural 
stability and high stress level implied that sharply reduced mechanical fatigue design life may be anticipated. 

These researchers have conducted a number of modeling studies of various stent types, geometries, and 
materials. These have included stent only, stent and catheter balloon only, and stent/catheter balloon/arterial 
wall. In general, good qualitative correlation has been achieved between the predicted deployed stent 
geometry obtained from the stent models and the observed deployed stent geometry obtained empirically. 
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Additional work is in progress to extend the connputational results for other stent geometries, types, materials, 
stent coating materials, effects of stent/balloon combinations, and evaluation of deployment strategies. 





This information is available in alternative formats upon request by individuals with disabiiities. Please send 
email to alt- format(g>msi, umn.edu or call 612-625-1818. 

f *Rt • htlp://www. msi. umn. edu:80/generat/BuUetinA/oL 1 4-No. 3/Juty98.html ^ 
This page last modified on Friday, 20'Aug-1999 09:17:18 CDT 
Website related questions or problems should be directed to webmaster<g>msiMmn.edu 
The Supercomputing Institute does not collect personal information on visitors to our website. For the University of Minnesota policy, 

see www.privacyMmn.edu. 



o 



4 



mechanics and 
engineering 

ELSEVIER Compui. Methods Appt. Mech. Engrg. 132 (1996) 45-61 ■■ 



XP-000972573 Oh 




Large strain analysis of soft biological membranes: 
Formulation and finite element analysis 

Gerhard A. Holzapfel^'*, Robert Eberlein*', Peter Wriggers^, Hans W. Weizsacker*" 

'Division of Applied Mechanics, Deparimeni of Mechanical Engineering, Stanford University, Stanford, CA 94M}S'4040, USA 
^Institui fur Mechanik, TH Darmstadt, Hochschuistrafie 1, D-63289 Darmstadt, Germany 
' Physiologisches Insiitut, Karl franzens Universiidi Graz, Horrachgasse 21, A-SOIS Croz, Austria 

Received 18 April 1995; revised 17 October 1995 



Abstract 

This paper presents a general formulation of thin incompressible n\embranes and investigates the behavior of soft biotissues 
using ihe finite eiemeni method. In particular the underlying hyperelaslic model is chosen to examine the highly non-linear 
constitutive relation of blood vessels which arc considered to be perfectly elastic, homogeneous and (nearly) incompressible. 
First, the stress-dcformaiion relation and the elastic tangent moduli are derived in a very general material setting which is 
subsequently specified for blood vessels in terms of Green-Lagrangian strains. Based on the principle of virtual work the finite 
element equations are provided and briefly discussed. Consistent linearization of the weak form of equilibrium and the external 
pressure term ensures a quadratic convergence rate of the iterative solution procedure. On the computational side of this work an 
effort was undertaken to show a novel approach on the investigation of soft tissue biomechanics. Representative numerical 
analyses of problems in vascular mechanics are discussed that show isochoric finite deformations (large rotations and large 
strains). In particular, a numerical simulation of the interaction between an inflated balloon catheter and a plaque deposit on the 
wall of a blood vessel is presented. 



1. Introduction and overview 

The objective of this study is to present a new numerical approach to the mechanical behavior of 
biotissues — in particular of blood vessels — which is based on the framework of non-linear continuum 
mechanics. One novel contribution to the field of computational mechanics is the Finite-Element (FE) 
simulation of the rai abdominal aorta and of the angioplasty, a technique in which an inflated balloon 
catheter interacts with a diseased blood vessel with the purpose to widen the lumen of the vessel. The 
numerical analyses have been managed on the basis of a hyperelastic biomaterial proposed by Fung and 
co-workers [1^3]. An implementation of the constitutive model within a finite element context has 
never been taken notice of in the literature, as far as the authors are aware. 

Often the attempt is made to simulate soft vascular tissues with isotropic thin-walled (latex) rubber 
tubes (see e.g. [4-6]), which is a reasonable first approximation to the actual material behavior if the 
blood pressure domain of arteries is considered far below the physiological range. In fact within the 
range of low pressures the artery reacts at small strains, displays Yubber-like' elastic behavior and may 
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be modeled, e.g. by means of the classical isotropic nco-Ho<'>kciui or Moonev-Rivlin pv)ienuals (see j7, 
8]). However, a blood vessel under physiological conditions and higher internal prcssiives behaves as a 
mechanical anisotropc (cylindricaborthoirope). primariix recognized b\ Patel and Frey which to 
model and lo analyze is our main goal. A rubbcr-bused constitutive model is certainly not an adequate 
material for simulating arteries over wide ranees of deformations. 

Vascular walls arc (1) (nearly) incompressible for in vivo strains (first shown hv Carew et al. [10|), 
(2) perfectly elastic (see e.g. [II]) and (3) show significant stiffening effects with increasing strains and 
internal pressures. In general, an artery is considered to be (4) homogeneous throughout the entire wall 
thickness and (5) free of residual stresses and strains, which in fact is disproved by the early work of 
Roach and Burton [12] (who indicate the significantly different mechanical behavior of elastin and 
collagen in the vascular wall), and by McDonald [13, p. 260] (who addresses an observation for the first 
time, namely, that an arterial strip which is cut through the thickness will open up), respeciivcly. 

Properties (1) and (2) are generally accepted in the literature, while the last two considerations are 
very open to question and still an area of ongoing research in physical modeling of biotissues. e.g. in 
recent studies the artery wall is regarded as a thrcc-conceniric-layered composite material (from. the 
lumen outward; intima, media and adventitia) with different mechanical properties in the zero-stress 
state (see [14, 15]). For recent studies on residual stresses (see e.g. [15-17]). however, the biomechani- 
eal importance of residual stresses and strains was already recognized [18]. 

Based on assumptions (})-(5) there are some phenomenological constitutive models in the literature 
with the same objective, namely, to correlate closely with experimental stress-deformation data for the 
entire occurring finite strain domain. The material models which are proposed in the form of more or 
less accurate strain (stored)-energy functions (potentials) are either given in a logarithmic form (see 
[19-21]), in a polynomial type (see [22-24]) or in the most commonly used exponential form presented 
in [14, 25-28] and in [1-3] which was also adopted in this paper. A selected review of the literature for 
constitutive laws of arterial walls was recently presented by Hayashi [21]. A review of a numerical 
approach in arterial wall mechanics with the FE method was given in a recent summarv paper bv Simon 
et al. (29). 

The theoretical side of the paper starts with a review of the basic equations for thin incompressible 
membrane surfaces. The formulation encompasses large strain elasticits and is exclusively based on the 
right Cauchy-Green tensor and the principal stretches as a strain measure. Explicit expressions for the 
stresses and the symmetric tangent moduli are directly derived from the stored-energy function and are 
provided in a form which is independent of any particular choice of a hyperelasiic material. Observe 
that for membranous structures it is sufficient to use only the tangential stresses in the deformed 
membrane surface (plane stress state). Within this two-dimensional approach the plane stress condition 
is imposed by introducing an additional constraint in the potential through an arbitrary scalar 
(Lagrangian multiplier) known as the hydrostatic pressure. The pressure which can always be expressed 
in closed form enforces the plane stress condition exactly. Next, we particularize the frame-invariant 
strain-energy function for analyzing the highly non-linear stress-strain response of an arterial wall. The 
consistent linearized tangent moduli is explicitly given as a novel result. 

With regard to a numerical simulation which is based on the FE method, the principle of virtual work 
is introduced. A consistent linearization of the weak form of equilibrium leads to the governing basic 
equations essential for a FE simulation. For problems involving highly deformable elastic membranes, 
the pressure load is generally deformation dependent, which must be accounted for in the formulation. 
Special emphasis is devoted to the linearization of the external pressure term. A spatial discretization of 
a continuous membrane represented by assembling a finite number of displacement-based finite 
elements leads to the tangent stiffness relation and to the global residual force (out-of-balance force) 
which ensures quadratic convergence within Newton's method. 

On the computational side of the paper, numerical results are illustrated for two selected axi- 
symmetrica! problems of soft-tissue engineering undergoing finite strains: the hyperelastic model was 
fitted to experimental data obtained from a rat abdominal aorta. The contact of a balloon catheter on 
atheromatous plaque deposits located at the innermost layer of the wall of an artery is particularly 
addressed. 
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2. Basic equations rff thin membranes with incompressible material behavior f 

In this section wc derive the basic cquiiiions used to describe ihe finiic strain behavior of thin 
membrane surfaces which are assumed to be hyporelastic and incompressible. For a more comprehen- 
sive studv on finite strain elasticity the reader may be referred to the standard books (e.g. [30-32|). 

The presentation of a hyperclastic constitutive relationship which approximates the mechanical 
behavior of blood vessels closes the section. 

2J. Kinematics of deformation and strain-energy function 

Consider a hyperclastic membrane (part of U^) undergoing large deformations (rotations and strains) 
which carry a material particle initially located at position A' = A'(6)") (t?" In -1.2 ^^^c curvilinear Gaussian 
coordinates) into a spatial position x = x{0") through the displacement vector k := a: - A'. Accordingly, 
we define a reference placement as ^„CiR'^ and a current placement (deformed configuration) 

as SSC(R\ 

The tangential vectors G„ and the external unit normal vector field A' (by the cross product) are 
represented as 

The Gaussian frame {G„.N} forms a right-handed vector triad. Let G "* be the reciprocal basis, 

such that 5f=(G.,.G^) and G^p = (G„,g/>. where Sf denotes the Kronccker delta, (•,-) the 
Euclidean inner product and G^^ the metric coefficients, respectively. Similar definitions are introduced 
for the current configuration and are subsequently denoted by lower case letters. 

The measure of deformation based on the two-dimensional theory is described by the symmetric right 
Cauchy-Green tensor C which is expressed by the material deformation gradient F. i.e. 

C:=fV=5„^G"®G^ , F :^Vx^g„^G^ . (2) 

According to the representation theorems for invariants (see Section 10 in [30]), the strain-energy 
function '7f'(per unit reference volume) for hyperelasiic materials may be formulated as functions of the 
strain invariants (see [30, p. 317]) 

v'=iiu^)=^nHi,j,j,). (3) 

satisfying the condition 7?X1> 1, i) = 0, where , 2..^ characterizes the three principal invariants of 

the right Cauchy-Green tensor C := C^^gG '* ^S^G*^, which are defined with respect to principal axes and 
a Cartesian basis as follows: 

=/,(C):-tr[G] = C^,+C, + C33, 

= /,(C) := y {(tr[C])' - tr[c']} = C^C. + CnQ:. + . W 

A-^A(C):=det[C] = C,,Q:C33. 
Equivalently, W can be formulated in the alternative form 



^/'■=7r(A^) = 7/U,,A2, A3) 



where {A"^} , 3 are the three real and positive (and non-zero) eigenvalues of the eigenvalue problem 
{C- X\\)e^ =o"and correspond to the three eigenvectors {^^}^ = ,.:.3 (principal directions of strains). 
The three values {A^},,^, , 3 are called principal stretches and 1 = G^ iS^G*'* denotes the identity tensor. 
The invariants of Eq. (4) expressed as functions of A^ are given as 
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V5 /, = a; + a; -f A] , /: = a;a_; + a;a? + a;a^ 



/, = a;a;a: 



(6) * 



Incompressible materials have the property thai the volume is preserved, which requires the 
following well-known condition (sec e.g. (30, p. 71]) 



or rewriiien in the form of principal stretches, with Eq. (6),, 



AiA,A, 1 



(J) 



where the dilatational ratio J characterizes the ratio between the incremental deformed volume di; of a 
body and its corresponding incremental undeformed volume dVA 

Through condition (S), the third principal stretch normal to the membrane degenerates to a 
dependent value, namely. 



1_ 

A, A, 



(9) 



Therefore, the current thickness of the membrane h is always determined by relation A. = h/ where 
// denotes the initial thickness. Now, with help of Eq. (9), the non-vanishing components of the right 
Cauchy-Green strain tensor C in terms of principal stretches are given in the convenient matrix setting: 



C,, 0 



0 

C„ 0 



a; 0 
0 a; 



0 
0 



(10) 



0 0 CjJ |_o 0 (\,\,y- = {h/H)- 

2.2. General fonmilaiion of the stress and elasticity tensor for incompressible membranes 

For incompressible hyperelastic continua the strain-energy density H' is given through the addition of 
a Lagrangian multipher p, i.e. 



(11) 



wherein 7r describes the isochoric deformation and p is an unknown scalar with the significance of a 
hydrostatic pressure, respectively. 

The stress-deformation relation for a Green elastic or hyperelastic material is generated from the 
continuous strain-energy function Eq. (11), by the relation 



dE ac ^ dC 



(12) 



where E = i{C - I) is the Green-Lagrangian (material) strain tensor and 5 : = 5^^C 0G ^ denotes 
the corresponding work conjugate second Piola-Kirchhoff stress tensor._ For a membrane, i.e. E: = 
f:,.pG"®G^ with definition 2E^^ : = - G,^ (see e.g. [33]) and S ^ dW/ dE := S^^C (S>G^ , respec- 
tively. With the partial (Frechet) derivative of quantity J relative to the svmmetric strain tensor C (see 
[30, p. 25]), i.e. 



^ dJ - - 1 
2—r = JC 

ac 



(13) 



and the incompressibility condition given by Eq. (7)^, relation (12), can be rewritten in the form 
(compare Fung [34, Eq. (12), p. 450]) 
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5 = 2 



dC 




(14) 



where C characterizes the inverse right Cauchy-Green tensor. 

With the use of the chain rule and Eqs. (9), (10), (14), a straightforward computation shows the 
principal values of the symmetric Piola-Kirchhoff stresses, which are given in the following matrix 
setting: 



1 dW p 
1_ 



1 



0 



0 



A, A, 



+ pA;A; 



• (15) 



Furthermore, using the plane stress constraint {S^y 
pressure to be explicitly expressed as 



P = ^ 



1 



A, A. aAjCAi, A.) 



0) the third row of Eq. (15) allows the hydrostatic 



(16) 



Observe that for membranous structures the incompressibility condition is already incorporated due 
to Eq. (9) and the Lagrangian multiplier p introduced in Eq. (11) serves here as an arbitrary scalar 
function to enforce the plane stress condition. 

Hence, a back-substitution of the pressure into relations (15),, (15). leads to the following two 
non-zero components 5'^lo-i,2» 



•^1 •"•^11 ~ 

5,:= 5,. = 



1 



1 dW 

A, aA, a-;a, aA,(A,, A2) ' 

1 dW 1 

A. A,A^ aA3(A,, A2) • 



(17) 



Through a rotational transformation, the principal stresses 5„ can be easily converted in stresses 
acting in any arbitrary position. Alternative stress measures are obtained by suitable transformations, 
e.g. the stress response in the Eulerian (spatial) description — the Cauchy-stress tensor a for a 
membrane — follows with a push forward operation and the incompressibility condition (7)^, i,e, 
o- -FSF"^, with F defined in (Z),- 

In the next part of this section the symmetric fourth-order tensor C — the elastic tangent moduli in the 
Lagrangian description— is derived by linearizing the second Piola-Kirchoff stress tensor. Differentiate 
Eq. (12)2 once more with respect to E and use the chain rule, then we obtain an explicit expression for 
the desired consistent incompressible' moduli: 



as 1 - 

dS = 2 — ::r-dC 
dC 2 



with 




(IS) 



in components, i.e. 
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Analogous lo the structure presented in Eq. (15). the components of the elastic moduli, i.e. (19). may 
be expressed in terms of principal stretches by means of Eqs. (9). (10), (13). Hence, with the plane 
stress condition, the elasticity tensor components for incompressible membrane shells are summarized 
in the following material matrix denoted by Z), i.e. 



D = 



A, dk, U, a; ■ 



1 fi/J 



A, i)k._ \A, dA, / 



1 bp 

aJaT^ 



J d_ /_]_ ^\ ] 

A^ rU, \A, f'A, / ^ AjA, 



dp_ 



A, f)A, ra J " A^ ^ a: cU, 



(20) 



Herein, /? is the hydrostatic pressure expressed in Eq. (16) and its derivatives with respect to A„|„^, > 
arc computed as 



1 cr>/' 



Ax 



where A, is a dependent value with regard to Eq. (9). 

The presented formulation is valid for all incompressible membranous structures which are based on 
hyperelastic materials and form an essential general basis for subsequent sections. 



2.3. Constitutive model problem of soft biological tissues 

With the derived concept presented in previous sections the continuum framework only depends on 
the specification of one scalar function, i.e. the strain-energy density the pursuit of which is the 
purpose of this section. 

In this approach the non-linear stress-strain behavior of biological tissues is described with the most 
concise potential used in literature, i.e. the exponential potential proposed by Fung et al. [1] and 
Chuong and Fung [2, 3] 



7/B,o = c{exp[G]-l} 



(22) 



where c denotes a constant with the unit of pressure. Within the concept of finite bio-elastomechanics Q 
denotes a quadratic function of the components of the Green-Lagrangian strain tensor. A particular 
form for blood vessels mioht be chosen as 



(23) 



where Oj, . . . ,0^ are non-dimensional constants characterizing the blood vessel and E^^ E^ and 
denote the three Green-Lagrangian strain components in the circumferential, longitudinal and radial 
direction of, e.g. an artery, respectively; E^g{A B) are assumed to be zero. 

The theoretical prediction of the hydrostatic pressure follows from Eq. (16) using Eqs. (22), (23) and 



relation 2E. 



2c 



P = - 



(A, A.) 



T [a^E^ + a^E, + a^E^] exp[Q] 



(24) 



The material coefficients o,,...,^^ in the stored-energy function can be identified from a fitting 
process with the goal to minimize the errors between known experimental data and theoretical 
prediction. 



3NSDOCID: <XP 972573A_I._> 



^ If a biopd vessel wall is treated as a membrane shell, only ihe average (uniform) strains in rhe 
circumferential (£,) and longitudinal (E.) direction arc of interest. These simplified assumptions lead lo 
the so-called two-dimensional approach, then, ihc analytical reprcseniaiion of nio- (--)• (23). for 

.a vessel scgmeni is obtained merely by setting and «„ to zero (see [I]). The alternative 

strain-energy function with four remaining material consianis basically represents a 'compressible' 
material model in two dimensions. 

The specification of the non-linear constitutive equation (17) and the consistent linearized tangent 
(20) results merely by computing the fust and second derivative of )/"ii,o with respect to the principal 
stretches, i.e. 

5, = 2r(«.£, + a,E,) cxpfQJj ^ [D., 

5, = 2c(a,£. + fl,£i)exp[<2lj ' L^2i ^-J ' 

with 

D,, =2c{fl, + 2(fl_,£, + fl,£,r} exp[Q) , (26) 
D,, = D„ = 2c{a, + 2(/7,£, + a,E,)(a.^E, + a,E^)} explQ) . 



3. Finite element formulation 

3.]. The pr'mcipie of viriitnl work 

With regard to a finite element concept the balance laws of a thin hyperelaslic membrane are 
enforced in weak form. Therefore, we apply the principle of virtual displacements, in which the 
displacement field u is taken to be a fundamental unknown on the membrane surface = ") 

of the boundary ^^^^ of the elastic continuum body. 

Let t '.^ p n be defined as a normal traction vector on the membrane surface f)5J;»„ with - 
c?^„„Ua^„^ and d^oi. ^ ^^0./ 0- '^^^ ^^^^ pressure p is a given constant per unit of the deformed 
differentia! surface element dj and n denotes the external unit normal vector field perpendicular to the 
loaded actual part Further, if 17 denotes a vector of the virtual displacements, then, the weak form 
of the balance equation may be formulated by the stalionarv value of the functional which is defined 
by 

^(«, Ti):=[ hE:SHdS-\ p<n, rj) ds = 0 , (27) 



for any admissible variations 17, which satisfying the mechanical boundary conditions tj]^^ = 0. In Eq. 
(27) E and S are the corresponding strain and stress tensors for a membrane as discussed in Section 2.2. 
The two terms of ^ identified as SWj^, and hV\^^ denote the internal and external virtual work, 
respectively. The second term, is the virtual work of the pressure on the boundary analyzed and 

extensively discussed by Schweizerhof and Ramm [35], Bufler [36] and with an application to 
axisymmeiric problems given by Simo et al. [37]. 

By using a relation similar to Eq. (1); for the current configuration and expressing the deformed 
differential surface element ds through the undeformed surface d5, the external virtual work of Eq. 
(27) may be rewritten to: 



(28) 
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^ 0\A. HolzQpfei ei at. I Comput. Methods AppL Mech. Engrg. J 32 (1996) 45-61 

The non-iincarity in u — caused both by the geometry and the material model — is implicitly present in 
the functional ^(w, tj). Therefore, with regard to an iterative algorithm of Newton's type, a consistent 
linearization of ^ at a known state u is necessary (see e.g. [38] or [31, Chs. 4-5]). 
The linearization of the functional is obtained by 

L<S{u, >), Au) := ^ + A'i^ + • • • = 0 , (29) 

where § denotes constant and linear terms, respectively. In this linearized problem the material 
linearization operator A{*} (Gateaux operator) is defined by the following rule (see (31, p. 185]): 

A(.}(«, 7,, A«) 0{-)(«, r?) • Au = ^ {•}{u + b A«)|,.o • (30) 

In view of the last two definitions (29) and (30) the functional ^{u, rj) together with Eq. (28) leads to: 

A^(ii, 7), Aii):=D^^(tf, tj)-Aiis ( (SE:A5-f5A£:5)//dS 



-I 



\G, X G,| ~ 17) d5 = - t)) . (31) 

This linearized problem must be solved for the increments An until the residual 17) on the 
right-hand side of Eq. (31) vanishes within a prescribed tolerance. 

3.2. Spatial discretization 

The domain 55 C K" is divided in an arbitrary number of isoparametric finite elements of the form 

"elm M 

^ = A,^^ where w^j^ is the number of elements of the region and A denotes the assembly 
operator, respectively. With each single finite element, ^'CS^'', the position vector Xy the displace- 
ment field u, the admissible variations ij and the incremental displacements A« are interpolated by the 
standard conventions of finite elements: 

"node* "node* 

(32) 

"nodei "nodes 

where the indicator h denotes the approximated function within the scope of an isoparametric concept. 
In Eq. (32), X„, v„, Av„ (« = 1, - ♦ ^ , /i^iodcs) represent the nodal values of A'\ u\ ri\ Au\ 
respectively. The number of nodes of a single element 58' is «^odej ^i"^ isoparametric shape 

functions associated with nodes n. 

In addition, we introduce the commonly used matrix containing derivatives of the shape functions. 
Then, the approximations for the virtual Green-Lagrangian strains SE and for the incremental strains 
AE yield: 

"nodes 

"V (33) 

"nodes 

With the above approximations, the discrete form of functional (29) is generated by 

"elm 

L<§\u, 7,, ^uf := A [^\v, q)\^^ + A^ Ai;)^ = 0 . (34) 
For one finite element, Eq. (31) becomes: 
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^ o./i. nOiZufjjct lu tu. I x^utnput. mvtnuui /utrt/i. c/ii;/^. tjj, (ijf^uy t.^-wi j_, 



-I 



»n. lev X Gjl 



^ ((A«: X - Agi X gll n') 65' 



= -[[ 5£*:5*//d5*-/ {gly^g'i.v'ydS"], (35) 

where SAfi'', AS'' and A^^ are approximate linearized functions of ihe virtual strain tensor, the second 
Pioia-Kirchoff stress tensor and the base vectors in the current configuration, respectively. 

For a FE formulation a compact matrix notation will be used subsequently. With the help of relations 
(32), (33) and the assembling operator A, the discrete functional given in Eq. (35) can be rewritten in a 
FE approximation as 

"^'^ /"Sodc* '"node, 

"elm /rt^oj^j \ 

= -A(S glK). (36) 

Herein, is the tangent stiffness matrix for an element and is determined by the initial displacement 
matrix B^DB^, the initial stress matrix K^^ and the crucial load tangent stiffness matrix caused by 
deformation dependent loads, as follows: 

f r ^ ^nodcs 

:= BlDB^H dS' + K„„ - J Ar„ S iN,.^N^.z - N,,,N^.,)i, dS' . (37) 

Ja&l •'a^o.r |0, XG2I p = l 

where N„ are the derivatives of the shape function N„ with respect to 0" and the known relation 

Ai;„xjc^=jc^Ai;^ VAr^GR" (38) 

between the cross product and a skew-symmetric matrix has been used. Observe the 2x2 material 
matrix /> of Eq. (37) which was explicitly presented by Eq. (20). 

On the right-hand side of Eq. (36), R'^ denotes the residual vector for the finite element which is 
defined by 

K ■-L,b:s'-H6S'- - d5^ , (39) 

where the first part denotes the stress divergence vector and the second part is the load vector including 
the deformation dependent loading terms. For convenience, in Eq. (39) we have introduced the 
definition 



::=g?x«5= 2 2 A/»a^„.2*,x*™. (40) 



"nodes ""nodes 

where x„ represents the nodal values of the current position vector x and contains here the discrete 
second Piola-Kirchoff stress components obtained via the hyperelastic constitutive law, i.e. = 
r'}^t-DBy^, 

Finally, for arbitrary virtual displacements ^ e R" a system of algebraic equilibrium equations follows 
from Eq. (36), which subsequently leads to the unknown corrector increments Ap. 
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4. Numerical examples 

Finally, we focus on two characierislic problems from the field of physiological and medical research. 
The objective of the present section is to demonstrate the performance of the preceding formulation, 
the phenomenological model and the numerical efficiency of the proposed finite element implementa- 
tion discussed above. A second main goal is to indicate a new practicable path for a numerical approach 
in the complex field of soft tissue mechanics. On the basis of a sound continuum formulation, the 
examples illustrate some novel numerical results on non-linear finite strain bio-elastomechanics using 
the finite element technique. 

Due to the axisymmelric geometry of the test problems with known principal axes of deformation, an 
axisymmetric two-node membrane element was employed which ensures quadratic convergence 
behavior within the Newton-Raphson algorithm. For the developed finite elements which account for 
large membrane strains it is sufficient to discreiize just the displacement field by linear interpolation 
functions (for more details on this subject see [37, 39]). 

The constitutive model, the finite element formulation and the numerical algorithm for deforrrration 
dependent loadings have been implemented in the finite element program FEAP, developed by R.L. 
Taylor and J.C. Simo (see [40] for a documentation). 

4.1. Inflation of a prestrained blood vessel — rai aorta 

The first example investigates the non-linear mechanical behavior of an abdominal aorta from a male 
Wistar rat via FE method. For the experimental set-up employed in this investigation and for a detailed 
discussion of the analytical assumptions and framework, the reader may be referred to [41-43] with 
more references therein. With regard to a numerical simulation it seems to be important to give a rough 
explanation of the experimental process. 

Basically, an aorta segment cannulated at one end and closed at the other was inserted between the 
clamps of a tensile testing machine and stretched to several levels of successive lengths. At each length 
the blood vessel, which is connected to a piston pump is inflated with a static pressure from 0 up to 
200 mm Hg (mercury) (26.66 kPa) acting along the luminal wall. For a certain longitudinal stretch ratio 

and an applied pressure p, measurements were made for the outer diameter d and the axial force. In 
addition, the outer diameter D and the wall-thickness H of the undeformed configuration of the vessel 
segment were measured. The experimental results are reported in [43). 

A straightforward computation with the incompressibility condition (9) (assumption that the arterial 
wall deforms iso-volumetrically) leads to the current wall thickness, i.e. h = H/A, D/d, Subsequently, 
within an analytical framework, which basically results from the equilibrium condition of thin pressure 
vessels, the non-linear stress-strain relationship caii easily be computed. The results for the physiologi- 
cal axial prestretch, which is approximately A^ = 1.60 are given in Figs. 1 and 2 and denoted as 
'experimental'. 

As a next step the experimental stress-strain curves were fitted with a non-linear least-squares 
algorithm over all experimental data points for A, -1.60. The material parameters obtained were 
c = 0.61163 10^ Pa, a, =0.68238, = 0*60558 and 0^ = 0.61163, with a^=^ a^^^ 0^, = 0, which goes 
along with the two-dimensional approach discussed in Section 2.3. The obtained coefficients minimize 
the sum of squares of the difference between known experimental data and predicted theoretical 
results. The fitted curves, labeled as *theoreticar are presented in Figs. 1 and 2. 

Based on these material coefficients the numerical simulations were performed with only one 
axisymmetrical membrane element and the specific geometry 0.122 mm (wall-thickness) and 1.648 mm 
(outer diameter) for the load free (undeformed) configuration was employed (see [43]). 

The first loading procedure, the prestrained state for p =^ 0 (A, = 1.60) was simulated followed by a 
second loading process which increases the pressure with incremental values of Ap = 10 mm Hg up to 
200 mm Hg, see Figs. 1 and 2. Therein, the radial displacements on the ends of the vessels are unlocked 
and for the inflation test the tube was fixed in a longitudinal direction. Note that the inner pressure is 
always normal to the surface, which significantly changes during the simulation. Since the stiffness 
changes due to the deformation dependent loads (see third term in Eq. (37)), large load steps can also 
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Circumrerentia] Green- Lagrange strain Ee 

Fig. 1. Inflation of a prcstrained blood vessel; circumferential second Piola-Kirchhoff sxtQ%s 5^ versus circumferential Green- 
Lagrange strain £^ at A. ^ 1.60, experimental data, theoretical and finite element solution. 



Material parameter: 




0 ' ' ' 1 1 ■ i 

-0.1067 0 0.2 0.4 0.6 

Circumferential Green-Lagrange strain Eb 

Fig. 2. Inflation of a prestrained blood vessel: axial second Piola-Kirchhoff stress S. versus circumferential Green -Lagrange 
strain at A. ~ 1.60, experimental data, theoretical and finite element solution. 



Table 1 



Inflation of a prestrained blood-vessel: norms of the residual for a typical increment load \p ~ 10 mm Hg — showing quadratic 
convergence behavior 



Iteration 


Step 1 


Step 2 




(90-» 100 mm Hg) 


(100^ no mm Hg) 


1 


4.49923 X 10"" 


4.58039 X 10'°* 


2 


3.79605 X 10"''* 


3.48725x10"*'' 


3 


2.17250x10"" 


1.65929X10*°' 


4 


7.22517x10*"' 


3.80461X10-'* 
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be considered and a quadratic convergence of the unbalanced forces near the solution point is always 
achieved — for two characteristic load steps see Tabic 1. 

REMARK 6.1. The non-linear stress-strain relationships (Figs, i and 2) show qualitatively good 
agreement with the theoretical values which motivates to discretize geometrically more complex 
problems (like bifurcations) found in the field of tissue engineering. This is a task for future work. 

However, as seen in Fig. 1 the experimental stress-strain relationship of a prestraincd abdominal 
aorta is typically 'S-shaped\ which cannot be replicated with the employed strain-energy function of 
exponential form. A potential which is capable to simulate the overall pressure domain more accurately 
was recently proposed by Weizsacker et al. [27, 28). This new potential for elastic and muscular arteries 
is composed of a neo-Hookean potential and an anisotropic potential as used in this paper. It allows a 
representation of the 'biphasic' behavior of arteries for the first time. An implementation of the 
improved structure of the constitutive model within the general concept presented in Section 2 appears 
to be straightforward, the goal of a foUowup paper. □ 

4,2. Contact of a balloon catheter with an atherosclerotic artery 

The last example is concerned with the mechanics of balloon angioplasty occurring in clinical 
applications. 

Percutaneous transluminal angioplasty (PTA) follows the primary role to enlarge the luminal 
(cross-sectional) area of a chronic atherosclerotic blood vessel by inducing and inflating a (latex) 
balloon dilatation catheter, leading to decreased resistance to blood flow. However, dilatation catheters 
generate increases of intramural stresses and strains, which also depend on the catheter design (see e.g. 
[44]), The most damage during balloon angioplasty occurs to the intima and media (see [45j), which can 
give rise to restenosis and embolic symptoms (luminal thrombosis) — overzealous dilatation can damage 
the contractile function of the smooth-muscle cells (see e.g. [46] among many others). One aim of the 
angioplasty technique is to produce controlled injury with a predictable healing response [47]. For 
further details regarding the general characteristics of PTA, the standard references by Becker et al. 
[48], Veith et al. [49], Chs. 21-22 are recommended. 

There has been only a few investigations into the mechanics of balloon angioplasty published in the 
literature (see e.g. [46, 50]), Analytical solutions of the non-linear mechanics of balloon angioplasty 
involving relatively complex geometry are hardly available, experiments are often difficult to perform. 
An efficient FE approach, which can handle the contact of a tiny balloon-tipped catheter with an artery 
narrowed by fatty deposits is one goal of this example — to the author*s knowledge still not available in 
the literature. Another purpose of the present study is to contribute to the knowledge in the mechanics 
of this complex process associated with finite strain elasticity by visualizing the global deformation 
behavior of the structure and the intramural stresses. 

Our attention is restricted to fully axisymmetrical problems of the model atherosclerotic vessel and 
the balloon catheter. Geometry and material parameter of the blood vessel — an abdominal aorta of a 
male Wistar rat — are adopted from the first example. The artery segment with a length of 2.5 mm was 
discretized by 50 axisymmetrical finite elements and prestrained up to 4.0 mm with a simulation process 
according to the first example, then, the current diameter for the middle surface becomes 1.353 mm and 
the wall-thickness is 0.086 mm. 

The centric plaque has a lumen with radius = 0.39 mm, as illustrated in Fig. 3. Furthermore, the 
assumption was made that the idealized plaque in the middle of the artery is circularly shaped with 
radius /?2 = 0.28mrh and the total length of the stenosis is 3.8 mm with a thickness on both ends of 
0.01 mm. The plaque geometry of an idealized model stenosis was discretized by 500 axisymmetrical 
finite deformation elements, which are based on a mixed FE formulation and a Hadamard material 
model. Basically, the compressible case of the neo-Hookean material was used, for further details see 
[51]. The mixed formulation circumvents the well-known locking response for the incompressible limit 
observed by a pure displacement formulation. The stress-free centric plaque is assumed to rest 
friction less on the artery wall. 
Little attention has been paid to the biomechanical properties of atherosclerotic tissues, which can 
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Fig. 3. Coniact of a balloon catheter with an atherosclerotic artery: discretization of an idealized atherosclerotic artery cross 
section with a balloon catheter showing the prestrained (reference) configuration. There are 650 axisymmctric elements: 50 artery, 
500 plaque and 100 catheter (a); close-up view of the discrelized centric model plaque (dimensions in miUimcier) (b). 



significantly vary from specimen to specimen, from hard to soft. Artery plaque systems behave as 
non-linear elastic imder the high stresses imposed by a percuianeously introduced balloon catheter. Very 
little experimental data regarding non-linearity* homogeneity and compressibility of the plaque within 
the physiological range are available (see e.g. (52J). The determination of material properties is still an 
area of incomplete knowledge and aKvays produces controversial opinions among scientists. For our 
simulations the assumption was made that the biomaterial of the atherosclerotic plaque is elastic and 
isotropic, described with typical value, i.e. the static shear modulus G — 500N/m" and the bulk 
modulus K = 724.6 N/m* (see [53, 54] for further details). Observe that the plaque generally is more 
rigid than the artery wall and most resistant to deformation. The strength of material composition of the 
dilatation balloon, which is held within a rigid catheter body as seen in Fig. 3, has to be sufficient to 
exert the pressure on areas of stenosis. For our simulations we have chosen a sphere with the diameter 
of the middle surface of 0.70 mm and a wall-thickness of 0.10 mm, which was discretized by 100 
axisymmetrical membrane elements. The catheter is modeled with the classical Mooney-Rivlin material 
(see [8]) which is adopted from empirical rubber-elasticity. The employed isotropic potential reads: 
'=lir ^ j^^^^ jyj^j^ „ 3) where the positive constants are chosen to be = 1/3.225 and /, denote the 

principal invariants of C as described in Eq. (4). 

Fig. 3 shows the system of the prestrained artery completely restrained at both ends with the cross 
section of the idealized plaque and the catheter, which is assumed to be the reference configuration. 
This is the point of departure for exertion of the dilating pressure on the stenosis. The atheromatous 
plaque is pressurized by the membrane surface of the balloon. The contact constraint is enforced by 
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means of a classical penalty method. The deformed shapes of the diseased vessel and the catheter with 
contour plots of radial Cauchy stresses in the plaque are displayed in Fig. 4 for a final state of the 
process (intraballoon pressure /? = -2.05 N/ mm", maximum contact pressure a, = -0.4 19 N / mm" ). 
The clastic catheter-plaque system which is pressed into coniaci by the applied pressure of the catheter 
give rise to a Hertzian contact pressure distribution, clearly seen in Fig. 5. 

Finally, the constitutive relations used have been implemented in a non-linear axisymmcirical 
quasi-Kirchhoff-type shell element, which is based on a formulation presented by Eberlcin'et al. 155). 
An analog contact simulation of the geometrical and physical non-linear behavior of the system were 
performed, which remarkably showed almost identical results. 

REMARK 6.2. It has to be mentioned that the final solution of this analysis will strongly depend upon 
the constitutive model used, which is probably the most poorly understood part in soft tissue mechanics 
and continually one goal of research. Implementations of efficient damage models are without doubt a 
next unalterable refinement towards a satisfactory constitutive modeling of vascular disorders. 

In future, 'computer aided balloon angioplasty* may help to predict indications for operations more 
precisely and may represent a method which enables different more efficient therapeutical interven- 
tions. Certainly the long-term goal of computer simulations of this important and very attractive clinical 
technique will be an aid to restrict occurrences of restenosis and to avoid difficulties by surgery which 
increases the rate of success for the treated patients. 

In this sense the present numerical study which is relevant to medicine should be considered as a 
basis for further needed research on the computational treatment of vascular surgery. □ 




Fig. 4. Contact of a balloon catheter with an atherosclerotic artery: deformed configuration of the balloon dilated atherosclerotic 
wall for the final state, showing Cauchy stresses in radial direction (dashed lines indicate the prestrained reference configuration). 
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Fij;. 5. ContncJ of a b;j|loon catheter with an nihemsclerotic artery: close-up view ol" the deformed system for the final state, 
showing the Hertzian contact pressure distribution in the pfaque (Cauchy stresses in radial direction). 



5. Conclusion 

In the present paper a classical constitutive model for biotissues has been employed to describe the 
highly non-linear stress response of homogeneous and incompressible hypcrelastic thin membranes — in 
particular of blood vessels. The condition of incompressibility has been imposed by a relation which 
expresses the radial deformation of the membrane through the in-plane (tangential) deformations. The 
assumption of a plane stress state has been enforced by introducing a potential through a Lagrangian 
multiplier (identified with pressure), which was explicitly determined. Derivatives of the strain-energy 
function with respect to Green-Lagrangian strains have led to the stresses and the elasticity tensor for 
an incompressible membrane. This framework is exclusively provided in a geometric setting relative to 
the reference configuration and encompasses a wide range of phenomenologically based hyperelastic 
constitutive models. To describe the perfectly elastic deformation behavior of soft biological (arterial) 
tissues, a particularization was carried out. 

On the basis of approximate displacement fields, a non-linear FE formulation for the finite strain 
behavior has been derived. A consistent linearization procedure of the governing equations allows the 
possibility to use large load steps by following quadraiically convergent equilibrium iterations. 
Consequently, the developed finite elements are very efficient and, in general, applicable for various 
materials in the large strain range. 

The performance of the described formulation and the accuracy of the numerical results have been 
demonstrated for two representative biomedical examples. For the first time in literature we were able 
to numerically simulate the mechanics of balloon angioplasty within the context of a FE method. The 
relevant information recorded in this study has been the connection between the applied internal 
pressure of the balloon catheter (loading) and the global deformation behavior of the structure (effect). 
The long-term trend in vascular surgery can be seen in accompanying computer-aided analyses which 
will help to visualize complex (patho-)physiological processes with the goal to minimize the risks for the 
treated patients. 
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Summary 

Endocoronary Ultrasonography and Angioplasty. 

G. Finet* At the end of the 20th century, endocoronary ultrasonography has become esta- 
blished in many catheter laboratories. The information provided by this technique of 
invasive imaging has many qualitative and quantitative features. Is endocoronary 
ultrasonography a research or an everyday clinical tool to be used systematically or 
occasionally in routine procedures ? 

Endocoronary ultrasonography is really a research tool because it is not redun- 
dant with conventional coronary angiography Ultrasonography investigates the 
whole arterial wall. Its spatial resolution is twice that of angiography. Its contribu- 
tion in this field is rich with many reports on arterial modelling, post-angioplasty 
remodelling, plaque rupture and atherothrombosis. the accurate description of the 
mechanical effects of angioplasty tools, two and three-dimensional quantification of 
atherosclerotic plaques and, finally, providing a very accurate biomechanical 
approach. 

Its use in routine procedures is only occasionally justified as there is no scientific 
proof in favour of a more widespread use with respect to changing the procedure or 
to evaluating patient prognosis. 

Endocoronary ultrasonography, a very rich technique, should be available in all 
catheter laboratories for occasional use. 

This technique is useful for ail cardiologists performing coronary angioplas^ by 
improving the evaluation of atheromatous plaques and countering the intrinsic 
ambiguities of coronary angioplasty. Arch Mai Coeur 1999 ; 92 : 1681-9. 

Resume 

A la veille de I'an 2000, Techographie endocoronaire est entree dans de nom- 
breux laboratoires de catheterisme cardiaque. Les informations fournies par cette 
technique d'imagerie invasive presentent de grandes qualites qualitatives et quan- 
titatives. Lechographie endocoronaire est-elle un outil de recherche ou un outll 
clinique a utiliser systematiquement ou occasionnellement en routine ? 

Elle est reellement un outil de recherche puisqu'il est non redondant avec I'an- 
giographie coronaire par rayons X, Lechographie explore toute la paroi arterielle. 
Sa resolution spatiale est 2 fois superieure a celle de I'angiographie. Sa contribu- 
tion dans ce domaine est riche avec de nombreux resultats sur le modelage arte- 
riel, le remodelage apres angioplastie. les ruptures de plaque et ratherothrom- 
bose. la description precise des effets mecariiques des outils utilises, les 
quantifications bidimensionnelles ou tridimensionnelies de la plaque d'atheroscle- 
rose, et enfin une approche blomecanique extremement precise. 

Son utilisation en routine ne peut etre qu'occasionnelle puisqu'il n'existe actuel- 
lement aucun element scientiflque suffisamment demonstratif pour autoriser une 
utilisation plus systematique en vue de changer soit nos procedures soit le pro- 
nostic de nos patients. 

Lechographie endocoronaire, technique d'une grande richesse, devrait etre pre- 
sente dans chaque laboratoire de catheterisme pour une utilisation occasionnelle 
en routine. 

Cette technique est utile pour chaque angioplasticien en enrichissant sa percep- 
tion de la plaque d'atherosclerose et en levant les ambiguTtes intrinseques de Tan- 
giographie coronaire. Arch Mai Coeur 1999 ; 92 : 1681-9. 
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sont actuellement disponibles : Tangiographle par 
rayons X et lechographie endovasculaire. L'imagerie 
des arteres coronaires par echographie endovasculaire 
est bien differente de l'imagerie obtenue par rayons X. 
II s'agit d une imagerie en temps reel, a haute resolu- 
tion, tomographique, de la lumiere et de la parol arte- 
rlelle. La question que I'on est en droit de se poser 
actuellement est la suivante : Fechographie endocoro- 
naire est-elle uniquement un outil de recherche ou un 
outil clinique de routine ? En d autres termes. son uti- 
lisation est-elle indispensable, pour la conduite dia- 
gnostique ou therapeutique de nos patients, en iabo- 
ratoire de catheterisme cardiaque ? 

Pour tenter de repondre a cette question, nous pre- 
sentons la methode d'utllisation actuelle de I'echogra- 
phie endocoronaire et nous tenterons d'apporter les 
reponses a ces questions au cours de la discussion. 

METHODES a . ...... . . 



Acquisition 

La tailie des catheters a ultrasons disponible a 
rheure actuelle va de 2.8 a 3 french. les frequences 
employees varient entre 30 et 40 Mhz. Les catheters a 
rotation m^canique sont plus utilises que les catheters 



electroniques. Effectlvement, le meilleur critere de 
qualite des catheters est represente par la surface du 
transducteur, environ 0.8 mm^. la puissance du fais- 
ceau ultrasonore emise est proportionnelle au carre de 
cette surface d'emission, Une miniaturisation exces- 
sive dlminue en coroilaire la puissance acoustique. 
principal determinant de la qualite des images echo- 
graphiques. Cependant. la qualite des catheters elec- 
troniques progresse. Les catheters actuels sont com- 
patibles 6F. 

routes les procedures d'exploration ^chogra- 
phique exigent finjection prealable de 250?mg d as- 
pirine, et d'un bolus de 2 000 unites d'hepirine. Lin- 
jection de dinltrate d'isosorbide ou de molsidomine 
doit etre realisee systematiquement avant la proce- 
dure. Le catheter a ultrasons est amene au-dela du 
segment arteriel que Ton desire explorer, sous 
controle radiologique, au niveau d*une zone precise 
et facilement reproductible telle qu'une bifurcation 
arterielle distale. Le retrait doit etre rendu automa- 
tique par Tutilisation d'un systemeautorise assurant 
la regularite du retrait et sa quantification. Une don- 
nee quantitative longitudinale est ainsi possible. 
Cette donnee intervient dans toutes les analyses 
volumiques quantitatives, dans la pertinence d'une 
reconstruction tridimensionnelle, et dans le repe- 
rage arteriel necessaire ^ la reproductibilite de la 






Fig. 1 - Differen^s aspects echographiques 
de la plaque d'otherosclerose coronaire. 
A : ospect echographiquement normal d'une 
artere coronaire avec Taspect lypique en 
3 couches caracterisant la limitanle elastique 
interne, soulignant la media, et caracterisont 
le port ad venticiel. B : plaque fibroceliulaire 
excentree. II est impossible sur des donnees 
echographiques de savoir si le contenu ane- 
chogene est lipidtque ou cellulaire (spu- 
meux, inflammototre, ou derive de produits 
sanguins]. C : plaque fibrocalcoire excen- 
tree. D : olaaue ulceree avec rupture cen- 
trole de fa plaque et comblement par un 
thrombus qui remplit lo niche formee par la 
rupture de plaque. 

Fjg. 1 - Different ultrasonic images of coro- 
nary atherosclerotic plaques. A: normal 
oppeo ranees or a coronary artery with three 
layers characterising the internal elostica 
underlying the media and the adventica. 
B: excentric fibrocellular plaque, tt is impos- 
sible to determine from an ultrasonic imaging 
if the anechogenic component is lipid or cellu- 
lar (spumous, inflammatory or blood pro- 
ducts). C: excentric fibro-calcific plaque. 
D: ulcerated plaque with central rupture filled 
with thrombus. 
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mesure d'un segment particulier, a 2 temps diffe- 
rents de Texploration. Toute ambigu'rte qualitative 
sur I'image du contour endoluminal, induite par la 
retrodiffusion acoustique du sang, doit etre levee 
par rinjection a fort debit d'un bolus de serum sale 
isotonique au travers du catheter-guide. 

Traitement qualitatif de Timage 

L'artere coronaire echographiquement normale 
est representee par un aspect caracteristique en 
3 couches, dent I'epaisseur est a peu pres 
constante. de 300 \xn\ environ. La detection d'une 
plaque d'atherosclerose est facile a identifier sur les 
images echographiques. La texture de I'image echo- 
graphiquG permet de distinguer 3 composants prin- 
cipaux (fig. 1). Une echogenicite moderee et homo- 
gene caracterise une structure d^ nature fibreuse, 
Une hyperechogenicite avec attenuation du faisceau 
ultrasonore traduit une zone calcaire. il peut y avoir 
quelquefois ambigui'te avec une fibrose dense de 
type hyaline. Une zone hypo- ou anechogene traduit 
une structure cellulaire. Elle n'est done pas syno- 
nyme de zone llpidique. Lhypo6chogenicite ou ane- 
chogenicit^ traduit Tabsence de reflecteur acous- 
tique. IJ s'agit done d'une structure riche en eau. Le 
contenu est done hpidique (extracellulaire) ou cellu- 
lalre (cellules spumeuses. cellules inflammatoires, 
produits sanguins). IMous avons decrit precedem- 
ment les arjibiguites diagnostiques propres a Tin- 
terpretation qualitative de Timage echographique 
endocoronaire [1]. 

Analyse quantitative 

Deux contours peuvent etre traces sur une image 
echographique endocoronaire : le contour endoluminal 
et le contour arteriel total. Le contour endoluminal 
consiste a suivre \e bord d'attaque de Techo engendre 
par interface sang-paroi arterielle. La detection de ce 
contour est le plus souvent realisable sans grande dif- 
ficulte mais quelques problemes specifiques peuvent 
apparaitre [2]. Le catheter a uitrasons est souvent en 
position excentree au contact de la paroi et peut gener 
la detection de cette surface qui disparatt du fait du 
contact avec le catheter. D^autre part rechogenicite du 
sang, surtout a 40 Mhz. peut etre telle que les impe- 
dances acoustiques du sang et de la parol deviennent 
tres proches. La visualisation de Tinterface disparait 
alors, genant la detection du contour. Une telle situa- 
tion apparait principalement lorsque le flux sanguin 
est ralenti, ce qui favorise la formation de rouleaux 
d*hematies, dont Timportance croit de maniere inver- 
sement proportionnelle a la Vitesse du sang, Ces rou- 
leaux deviennent alors des reflecteurs acoustiques 
Importants et genent ainsi considerablement le tra- 
gage des contours endoluminaux. Une methode simple 
permet d'ameliorer la detection visuelle de ce contour 
endoluminal et consiste a enlever les reflecteurs san- 
guins en injectant au travers du catheter porteur un 
bolus de serum sale isotonique a fort debit. Uutilisa- 
tion de produit de contraste acoustique est aussi une 
solution [3]. 
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Le second contour est le contour arteriel total et cor- 
respond au tragage effectue en dedans du bord d'at- 
taque de i'echo produit par finterface media, limitante 
elastique externe et adventice. Le trag age de ce contour 
est generalement facile mais peut dans certains cas etre 
partieliement absent. Cette disparition localisee du 
contour est due aux zones d 'absence de transmission 
des faisceaux ultrasonores par des composants hyper- 
reflecteurs de la plaque. Ces composants sont les 
plaques fibreuses denses et les plaques calcaires. 11 
convient alors d'extrapoler ce contour en fonction du 
contour visible en s aidant des images adjacentes a celle 
sur laquGlle on effectue les mesures. Au total 2 traces 
sont realises et fournissent 2 contours. Ces 2 contours 
creent 2 surfaces : la surface endoluminale et la surface 
arterielle totale. La difference de ces 2 surfaces fournit la 
surface intima + media ou plaque + media. Au total. 
3 surfaces sont ainsi obtenues par 2 tragages de 
contour. II Importe de detecter les deformations geome- 
triques engendrees par les changements de position du 
catheter uitrasons au sein de l'artere exploree, Len- 
semble de ces traitements quantitatafs doit etre dispo- 
nible sur une station de travail informatique dediee. 
conviviale et interactive (logidel IVUS 153.2, I6DP, Paris). 
La reconstruction tridimensionnelle et la visualisation 
(fig. 2) offrent un potentiel majeur dans Tanalyse quanti- 
tative volumique [4]. Lanalyse des artefacts de forma- 
tion de cette Imagerie tridimensionnelle est realisee afin 
d'assurer la validite des informations [5]. 




Fig, 2 ' Reconsfruction fridimensionnelle d'une stenose serree de I'in- 
terventriculaire anterieure proximale avec visualisation d'une 
branche loteroie (premiere septole). La texture de Timage echogro- 
phiaue est conservee sur ce type de reconstruction permettant une 
analyse longitudtnale de cette plaque fibro^cellulaire. 
Fig. 2 - Thr^-dimensional reconstruction of a severe proximal LAD 
stenosis with visualisation of a lateral branch (first septal). The texture 
of the ultrasonic image is preserved in this type of reconstruction, 
enabling longitudinal analysis of this flbrocellular plaque. 
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Pieges dans I'interpretation des images 

L'image echographique peut presenter des arte- 
facts, autrement dit des aspects specifiques, pro- 
duits par ie systeme d'imagerie qui ne refletent pas 
la realite [6. 7]. Ces pieges sont surtout induits par 
des artefacts formes sur Timage echographique 
endovasculaire iors de I'exploration des stents. Les 
stents, structures metalliques, sont hyperechoge- 
niques. Leur visualisation sur Pimage est done parti- 
cuJierement bien visible. C'est parce qu'eJles sont 



bien Vi^ibJes que ces structures ont ete largement 
exploitees et ont rendu semble-t-il I'echographie 
endovasculaire indispensable iors du depioiement de 
ces protheses. L'echographie endovasculaire semble 
ainsi benefique en evitant les implantations insuffi- 
santes, les defauts d'apposition des filaments sur la 
paroi arterielle. et enfin en fournissant des donnees 
quantitatives afin d optimiser ia taille du bailon et 
des pressions d^inflation [8-10]. Cependant, plusieurs 
points critiques apparaissent. Depuis Tintroduction 
de Techographie endovasculaire, Ja plupart des 




Fig. 3 - Demonstration des artefacts caracte* 
risttques de rexplorotion echographique 
endovasculaire des endoprotheses coro- 
naires. A gauche, les images sont acqui&es 
in vitro, sur Ie meme stent (Polmaz Schatz - 
PS 154) deploye a un diometre tres precis 
de 4 mm. A drotte, ('image est issue d'une 
simulation numerique complete d'un stent 
virtuel deploye d un.^iametre de 4 mm. 
A : le Iransducteur ultrasonore est centre, les 
contours sont parfaitement circulaires, rin- 
terpretation est optimale, les quantifications 
sont justes. B : le catheter est excentre au 
contact du stent mais I'axe longitudinal du 
cotheter est parallele a I'axe lonqiludinal de 
I'artere. Le stent n'est que peu deforme glo- 
bolement, cependant. les echos des fila- 
ments deviennent obliques pour certains 
d'entre eux et la largeur de ces echos varie 
en fonction de la distance aui les separe du 
Iransducteur. L'analyse qualitative p«ut etre 
ainsi perturbee et les quantifications sont 
moderement alterees. C : le transductgur est 
incline de 30 * a rinferieur du stent. Une 
deformation geometric|ue lypique en ellipse 
apparoTt. S'y associe les variations de for- 
geur des echos des filaments et leur change- 
ment d'orientation. L'ensemble de ces modi- 
fications sur l'image altere sigmficativement 
I'interpretation que t'on peut avoir sur {'im- 
plantation du stent et modifie profondemenf 
les quantifications. 

Fig. 3 • Demonstration of characteristic arte- 
facts of endovascutor ultrosonic imaging of 
coronary stents. On the left, the appearances 
acquired in vitro of the some stent (Palmaz 
Schatz - PS 154) expanded to a very precise 
diameter of 4 mm. On the right, on imoge 
from total numeric stimulation of o virtual 
stent expanded to 4 mm. A: the ultrasonic 
transducer is centred, the contours are per- 
fectly circular, the interpretation is optima) 
and quontification is occurote. B; the catheter 
is ex- centred in contact with the stent but the 
lonqitudincl axis of the catheter is parallel to 
the longitudinal axis of the artery. The stent is 
only slightly deformed globally ; however, the 
echoes of the filaments become oblique ond 
the width of these echoes varies occording to 
the distance between them and the transdu- 
cer. Qualitative analysis can be affected ond 

Quantification may be approximate. C: trans- 
ucer inclination at 30° within the stent. Typi- 
cal ellipsoid deformation. Associated varia- 
tions in the width of the echoes of the 
I'ilamenls and change in orientation. The 
changes significantly alter the interpretation 
of the stent implantation and greotly affect the 
quantification. 



1684 



ARCHIVES DES MALADIES DU CCEUR ET DES VAISSEAUX, lome 92. n' 1 1 (suppt), novembre }999 



MSDOCID: <XP 1014246A l_> 



^^..^ ECHOGRAPHIE ENDOCORONAIRE ET ANGIOPLASTIE 



Stents depJoyes (85 % environ) paraissent mal 
impiantes. Ces mauvaises implantations sennblent 
fondees sur plusieurs elements propres aux images 
echographiques : une geometrie asymetrique du 
stent a Torigine de sections elliptiques [8. 9], des fila- 
ments obliques, dont le nombre caracteriserait une 
expansion asymetrique du stent [9]. une geometrie 
parfois complexe en forme de petale ou les filaments 
semblent converger vers ie transducteur [8, 11], et 
enfin des defauts d'apposition de certains filaments 
centre la parol arterieJIe [8, 9]. Malheureusement. 
plusieurs paradoxes apparaissent dans la litterature 
et ne sont pourtant pas mis en exergue. Pourquol 
tous les defauts d'implantation ne peuvent etre cor- 
riges ou prevenus par une Implantation guidee par 
I*echographie endovasculaire ? Pourquoi malgre Tuti- 
llsation de hautes pressions d*inflation ou de ballons 
de diametre superieur, les criteres echographiques 
initialement retenus en faveur d'un defaut de 
deploiement sont ils toujours presents ? 

Nous avons analyse les interactions complexes qui 
lient les filaments metalliques des stents et leurs 
dispositions geometriques tridimensionnelies avec 
les proprietes acoustiques du catheter a ultrasons et 
sa position en realisant des etudes experimentales 
in vitro et des simulations numeriques completes 
(fig. 3) [12, 13], Lechographie endovasculaire appa- 
rait globalement comme une technique d*imagerie 
relativement inadaptee a Texpioration serieuse des 
stents. Effectivement. finterpretation des images 
echographiques omet completement les principes 
de formation'de cette image : position intraluminale 
ultrasonore, rotation du faisceau ultrasonore et 
interaction acoustique entre ies filamepts et le 
transducteur. Tous ces elements sont a Torigine 
d'artefacts qui perturbent nos interpretations quali- 
tatives et alterent significativement nos quantifica- 
tions. II apparait ainsi logique d'etre confronte a des 
discordances etonnantes entre les donnees echogra- 
phiques et angiographiques. 



DISCUSSION...^..--. - .. .... 

Une discussion sur Techographie endovasculaire doit 
s'articuler autour de 2 questions : Techographie endo- 
coronaire est-elle un outil de recherche ou est-elle un 
outil clinique a realiser au quotidien [1 4]. 

Outil de recherche ? 

L'echographie est certainement un outil de 
recherche de grande quality pulsqu'il est essentielle- 
ment non redondant avec la seule imagerie que Ton 
peut utiiiser pour J'exploration des arteres coronaires : 
Tangiographie par rayons X. Uechographie explore 
toute la paroi arterielle ce que ne fait pas Tangiogra- 
phie. E!le fournit des donnees quantitatives 2 fois plus 
precises que Tangiographie (90 [im versus 200 fim) 
[7, 15]. Ces 2 proprietes rendent cet examen indispen- 
sable pour Texploration in vivo de Tartere normale et 
surtout de la plaque d'atherosclerose. D'autre part, 
Techographie associe I'analyse morphologique de la 



plaque (bidimensionneJIe ou tridimensionnelle) et 
Tanalyse relativement precise de sa composition : ces 
2 caracteristiques sont observees en condition de 
charge physiologique. L'apport est considerable sur 
nos connaissances de revolution de la plaque d'athero- 
sclerose, sur le modelage et le remodeJage post-angio- 
plastie de cette plaque. Nous avons acquis la 
connaissance de 2 modelages ou remodelages carac- 
teristiques : la dilatation arterielle dite compensatrice 
decrite initialement par Glagov et coll. [16] et la 
constriction arterieJIe chronique. Cet aspect bimodal 
du modelage arteriel a ete decrit tant apres angioplas- 
tie que sur des arteres coronaires natives [17. 18]. 
Nous avons pu montrer que ces 2 modelages ne sont 
pas exclusifs et coexistent tres souvent (fig. 4. page 
suivante) [12]. De plus, le phenomene de constriction 
arterielle. initialement retrouve dans le remodelage 
post-angioplastie dans 25 % environ des cas est pre- 
sent de maniere focale ou giobale dans plus de 75 % 
des stenoses coronaires (fig. 5. page suivante). Les 
connaissances nouvelles ont fortement contribue a 
modifier notre perception de la plaque d'atheroscle- 
rose et ont eu des impacts importants : comprehen- 
sion de I'existence d'angiogrammes normaux alors que 
Tath^rosclerose est deja diffuse, critique sur la quanti- 
fication angiographique de Tatherosclerose, critique 
sur les segments de reference lors de la quantification 
des stenoses, amelioration dans la comprehension de 
Taction des outils d'angioplastie. Uechographie peut 
nous fournir dans certaines occasions anatomiques 
des informations sur Tadventice. tunique fort impor- 
tante dans les modelages pr^cedemment decrits, et 
qui subit elie aussi un modelage pathologique. Nous 
avons pu montrer que Tadventice en rapport avec 
['existence d'une plaque peut avoir une epaisseur jus- 
qu*a 4 fois superieure a celle d'une adventice d'artere 
normale [19], Uechographie s'est revelee beaucoup 
plus sensible dans la detection des calcifications [20] 
et plus precise dans la determination et la quantifica- 
tion des lesions coronaires. Elle s'est averee done une 
methode de choix pour revaluation des stenoses 
avant et surtout apres angioplastie. De plus, 
plusieurs etudes ont demontre que fechographie 
endovasculaire permettait une detection precoce de 
fatherosclerose acceleree du greffon en visualisant un 
epaississement intimal des vaisseaux coronaires angio- 
graphiquement normaux chez les patients explores a 
intervailes reguliers apres une transplantation car- 
diaque. Uanalyse relativement fine de la composition 
de la plaque a permis de montrer que dans Tangor 
instable, les patients avaient moins de plaques hyper- 
echogenes fibrocalcaires (10 %) et beaucoup plus de 
plaques a dominance cellulaire {90 %) done potentiel- 
lement riches en lipides et en cellules inflammatoires. 
La detection des ruptures spontanees ou de Tathero- 
thrombose est un avantage certain de I'imagerie echo- 
graphique endovasculaire. 

Uechographie comme outil de recherche permet 
d'aller plus loin dans I'extraction des donnees avec 
en particulier une analyse des comportements bio- 
mecaniques de la plaque d'atherosclerose directe- 
ment visualises chez le patient, Cela est le fruit 
d'un travail en collaboration avec des biomecani- 
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Fig. 4 - Exemple de remodelage constricHf 
arteriel. A : angiogramme d'une coronaire 
droite ou Ton peut voir un aspect ampullaire 
sur le deuxieme segment de cette arfrere. 
B : coupe echographique sur le segment 
proximal adjacent a Kampoule : pKenomene 
de modelage constrictif chronique majeur 
ovec reduction de 40 % de la taille totale de 
I'artere. C : coupe echographique endocoro- 
naire au niveau de I'ampoule. Aspect sub- 
normal avec epaississement fibreux intimal 
excentre. La masse atfierosclereuse est iden- 
tique dans ces 2 cos. Le processus stenosont 
global des segments ad|acents a rompoule 
n'est done induit que par un processus 
constrictif et non par un epaississement 
atherosclereux. 

Fig. 4 - Examples of constrlcMve arterial 
remodelling. A: Angiogramme of a right 
coronary artery with the ompullary oppea- 
ranees of the second segment. B: ultrasonic 
section oF the segmenf proximal to the 
ampulla : phenomenon of chronic constrictive 
modelling with 40% reduction of the total 
lumen of the artery. C: Endocoronary ultraso- 
nography at the ampulla. Subnormal appea- 
rances with eccentric fibrous intimal tnicke- 
ning. The atherosclerotic mass is the same in 
the two coses. The stenotic process of the 
adjacent segments to the ampulla is only 
induced hy a constrictive process, not by 
athrosclerotic thickening. 



ciens. Ces travaux biomecaniques nous apportent 
des elements riches d 'informations sur la. vuln6rabi- 
lite de la plaque d'une part et sur les forces com- 
pressives que Tartere fait subir a une eventuelle 



endoprothese (fig. 6). Enfin, Tanalyse volumique, 
permise grace a I'^chographie endovasculalre. per- 
met des analyses tres fines de la plaque d'athero- 
sclerose ou de la proliferation n6o-intimale intras- 




Fig. 5 - Aspect bimodol des modetages de 7 stenoses coronaires atherosclereuses natives. A : dHalation compensatrice defcillante puisque le proces* 
sus a malgre tout abouti a la formation d'une stenose endolumincle. B : constriction arterielle chronique globale particulierement bien objectivee. 
Ces 2 lesions sont des lesions sfenosantes serrees de 2 arteres interventriculaires onterieures, II s'ogit de reconstruction tridimensionnelte in vnro. 

Fig. 5 - Bimodal appearance of modelling of two native atherosclerotic coronary stenoses. A: inadequate compensatory dilatation as the process has 
resulted in the formation of endoluminol stenosis. B: chronic global arterial constriction, particularly well demon strotecf. These two lesions ore severe 
stenoses of the left anterior descending arteries. This is an in-vivo three-dimensional reconstruction. 
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tent (fig. 7, page suivante). Plusieurs grands tra- 
vaux sur la progression ou regression de la plaque 
d*atherosclerose sous statine vont connmencer en 
utilisant la quantification volumique apres recons- 
truction tridinnensionneile de la plaque d'athero- 
sclerose sur un long segment arteriel a des temps 
d'evolution differents. 

Uechographie endovasculaire est done un remar- 
quabJe outil de recherche, non redondant, precis, 
reproductible, et dont la contribution durant ces 
5 dernieres annees a ete reellement considerable. 

Imagerie clinique de routine ? 

La question que Ton doit se poser pour la pratique 
quotidienne de Techographie est la suivante : Techo- 
graphie doit-elle etre realisee systematiquement en 
cours d'angioplastie ou occasionnellement. En d^autres 
termes, existe-t-il actuellement des elements scienti- 
fiques qui nous obligent a realiser cette technique 
pour modifier directement, et pour chaque patient, 
une therapeutique ou un eventuel pronostic [21]. 

Uechographie est-elle necessaire pour selection- 
ner les outils d'angioplastie ? De nombreuses etudes 
largement dominees par les travaux du Washington 
Hospital Center, preconisent Tutilisation systema- 
tique de Kechographie avant angioplastie pour choi- 
sir VouiW adapte a la lesion traitee [22]. Les resul- 
tats cliniques qui en decoulent ne sont guere 
probants. ils ne modifient ni le resultat ni le pronos- 
tic des patients, lis sont a I'origine d'une augmenta- 
tion du temps de procedure, ils sont souvent lies a 
rutilisation de catheters a 8 a 10 french pour per- 
mettre rutilisation de tous les outils disponibles. 
Que dire d'une methode aussi lourde alors qu'en 
France par exemple, nous realisons nos angioplas- 
ties en 6F. avec de tres faibles doses d'heparine 
(2 500 UI). ballon et stent (dans 80 % des proce- 
dures), une ablation immediate de rintroducteur, et 
le tout avec une temps de procedure de 25 min 
environ. Le taux de thrombose aigu§ ou subaigue 
est inferieur a 1 .5%. 

L'echographie est-elle predictive de la restenose ? 
Les etudes sont tres divergentes et ambigues [22-25]. 




Fig. 6 - Analyse in vivo biomecanique du comportement de la plaque 
d'a^herosclerqse grace a I'associaHon de I'echogrophie endocoro- 
noire et de V<inafyse par les elements finis. A : image d'une plaque 
fibro-calcaire.avanf angioplastie. B : carrographie des conh-aintes de 
cisaillement ou apparait une importante lone de contrainte maxi- 
male partant de la lumiere e\ se prolongeanr vers la plaque calcaire. 
il s'agit d'une zone potentielle de rupture si les contramtes appliquees 
a I'artere deviennent superieures aux contraintes de rupture des ele- 
ments fibreux constituent la plaque. C : I'image echogrophique apres 
angioplastie objective une rupture de plaque. Cette rupture etait par- 
faitement previsible sur les iriiages de cartogrophie biomecanique. 
Fig. 6 - In-vivo biomechonicai analysis of the otheroscleroHc plaque by 
the association of endocoronory ultrasonography and analysis of the 
elements. ' ' i. . - 




mopping 

extending 1.^ ^ — , - — 1~' s ' t 

tiol rupture if the stresses on the artery become qreoter than the stress of 
rupture of the fibrous elements of the plaque. C: ultrasonographic ima- 
ging after angioplasty show^ing plaque rupture. This rupture wbs^per- 
fecHy predictable from the biomechonicai mapping data. 
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Fig. 7 - Exemples d'onalyses quanli^atives 
volumiques par echoprapnie endocoronaire 
d'une proiiferaHon neo-intimalega I'interieur 
d'un stent tubulaire. A : une vfie echogra- 
phique bidimensionnelle de la proliferation 
neo'intimale. B : reconstruction volumique 
tridimenstonnelle du volume intrastent. 
C : volume apres reconstruction tridimen- 
sionnelle de la proliferation neo-intlmale a 
I'interieur du stent. D : d'apres decoupage 
electrontque longitudinal de cette prolrfero- 
tion, il est possible d'en apprecier Tunifor- 
mite de la repartition au sem de ce stent. 
Fig. 7 - Examples of quantitaUve voiumic ana- 
lyses of endocoronary ultrasonography of 
neo-intimal proliferation with a tabular stent. 
A: two-dimensional ultrasonic appearances 
of neo-intimal proliferation. B: Three-dimen- 
sional voiumic reconstruction of intra-stent 
volume. C: volume after three-dimensional 
reconstruction of neo-intimal proliferation 
within the stent, D: after longitudinal electro- 
nic section of this proliferation, the uniformity 
of the distribution of this prolfferation within 
the stent may be observed. 



^importance de la plaque residuelle apres angioplastie 
et les importantes dissections semblent etre des fac- 
teurs predictifs de restenose. D'autres etudes ne 
retrouvent aucun critere predictif significatif [26. 27]. 
SI nous etions capables de predire une eventuelle 
restenose. le probleme qui se poserait a nous serait de 
savoir qu'elle serait alors notre action therapeutlque 
specifique. II n'y en a pas done cette prediction n'a 
guere d'interet en routine clinique. 

L'echographie est-elle un Indicateur precis de bonne 
implantation des stents? Ce grand sujet apparu en 
1994 repose sur I'interet potentiel de I'echographie 
endocoronaire pour Toptimisation du deploiement des 
endoprotheses fondee a partir de criteres issus des 
images echographiques reafis^es lors de I'implantation 
[28]. Certains auteurs ont observe que dans plus de 
80 % des cas, les stents etaient mal implantes [8]. II y 
a done discordance apparente entre I'aspect angiogra- 
phique et I'aspect echographique [29]. Ceia est un pre- 
mier paradoxe. La reaction immediate a cette observa- 
tion a ete d'utiiiser des ballons de taille superieure ou 
surtout de plus haute pression d'inflation. Les hautes 
inflations sont nees. Alors. un deuxieme paradoxe est 
apparu. Ce paradoxe est decrit dans quasiment toutes 
les publications importantes et souligne le fait que les 
criteres qui ont servi a la description d*une mauvaise 
implantation sont systematiquement retrouves apres 
I'utilisation de hautes pressions d'inflation. Nous avons 
pu expliquer puis demontrer par simulation numerique 



que rhyperechogenicite des filaments est a Toriglhe de 
nombreux artefacts qui modlfient tres sensiblement 
notre interpretation et aussi les quantifications. Ce que 
nous voyons sur I'image n'est pas le stent ou les fila- 
ments mais les echos du stent ou des filaments. La 
forme de ces echos est tres particuliere et les deforma- 
tions geometriques induites par les deplacements du 
catheter sont a I'origine de toutes ces deformations. II 
est done bien normal de retrouver ces memes criteres 
quelle que soit la surinflation imposee au stent. Quatre- 
N^ngt-quinze pour cent des etudes presentees utilisent 
toujours le meme stent : le Palmaz-Schatz. Nous avons 
pu montrer aussi par simulation numerique que la 
pression d'inflatlon necessaire pour ce stent est d'envi- 
ron 1 4 bars, elle est done en accord avec la litt^rature 
en general [30]. Mais cette pression n'est pas du tout 
applicable a I'ensemble des stents dont les proprietes 
mecaniques sont bien differentes. II n'y a pas de pres- 
sion magique mais une pression d'inflation a adapter 
en fonction du ballon utilise (pression nominale du bal- 
lon), de la pression de deploiement du stent et de son 
recoil elastique intrinseque (proprietes mecaniques du 
stent) qui doit etre compense par une surinflation 
determinee d'apres les courbes compliance fournies 
par le constructeur. Enfin. Schiele et coll.. au moyen 
d'une belle etude randomisee comparant 2 groupes 
d'angioplastie avec implantation de stents guides ou 
non guides par echographie endocoronaire, ne trou- 
vent aucune difference significative en terme de reste- 
nose [31]. 
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CONCLUSION . - ' 

L'echographie endocoronaire est une technique d'ima- 
gerie invasive qui nous fournit des informations d*une 
grande precision, completes, et non redondantes si on 
les compare a I'angiographie coronaire. C'est avant tout 
un outil de recherche dont les resuitats ont fortement 
contribue a modifier notre perception de la plaque 
d'atherosclerose. L'echographie ne peut etre a Theure 
actueiie un outil d'imagerie realise au quotidien et syste- 
matiquement au cours des procedures intervention- 
nelles. II est couteux et surtout aucun resultat probant 
ne nous conforte dans Hdee qu'il faille appliquer cette 
technique a tous nos patients au cours des angioplas- 
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ties. Cependant, cet outil d'une tres grande hchesse 
devrait devenir le « compagnon » indispensable de tout 
angioplasticien et devrait done etre present dans tous 
les laboratoires de catheterisme. Effectlvement Fecho- 
graphie endovasculaire pourrait etre compares a « Toeil 
de ranatomopathologiste » tant il est vrai que les infor- 
mations apportees sont comparables a celles du patho- 
logiste dans leur precision et dans leur justesse d'obser- 
vation, Uangiographie coronaire est d'interpretation 
quelquefois difficile et parfois ambigue : il est alors 
necessaire d'utiliser cette technique complementalre 
pour redresser notre interpretation angiographique tout 
en ameliorant nos interpretations angiographiques avec 
un gold standard anatomique. 



MOTS CLES : echogrophie endovasculaire, stenl, angioplastie, remodelage arteriel, atherosclerose. 
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